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INTRODUCTION. 
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A compelling  need  for  knowledge  of  the  properties  (temperature,  density, 
velocity,  composition,  etc.)  of  the  flows  produced  by  propulsion  systems 
has  recently  surfaced  in  many  areas.  The  need  has  arisen  because  of  the 
increased  effort  to  obtain  higher  energy  conversion  efficiencies  from 
combustion  processes,  to  define  visibility  and  other  observable  properties 
of  propulsion  system  exhausts,  and  to  make  quantitative  assessments  of 
environmental  effects  caused  by  exhaust  effluents.  One  area  of  concern 
is  with  rocket  propulsion  systems.  Of  specific  interest  to  the  Air  Force 
Rocket  Propulsion  Laboratory  (AFRPL)  is  knowledge  of  properties  within 
combustion  chambers  (particularly  for  solid  fueled  motors  employing 
metallized  propellants)  and  properties  affecting  exhaust  plume  observables 
(e.g.,  smoke,  infrared  radiation  signatures,  ultraviolet  radiation  sig- 
natures, radar  cross  sections).  Methods  of  making  measurements  within 
combustion  flows  must  avoid  influencing  the  very  phenomena  to  be  measured 
and,  generally,  must  be  made  from  without  (nonintrusive)  because  of  the 
harsh  environment.  The  field  of  measurement  associated  with  high  temperature 
gas  and  gas-particle  flows  is  generally  referred  to  as  "diagnostics."  This 
report  is  aimed  at  defining  the  needs  for  "combustion  diagnostics"  associated 
with  rocket  propulsion,  assessing  the  status  and  utility  of  the  applicable 
techniques,  and  establishing  a basis  for  the  research  and  development 
required  to  provide  techniques  to  meet  the  defined  needs. 

A few  definitive  statements  need  to  be  made  to  isolate  the  class  of 
measurement  problems  to  be  addressed  in  this  document.  Measurements  may 
be  categorized  in  the  following  ways: 

(1)  Ruler  - The  dimension  of  the  variable  to  be  measured 
and  the  observable  are  the  same,  usually  with  a one- 
to-one  correspondence  (e.g.,  ruler,  volume  vessel, 
protractor) . 
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(2) 


Gauge,  meter  - The  dimension  of  the  variable  to  be  measured 
is  transformed  to  a mechanical  indicator  or  dial  for  obser- 
vation (e.g.,  pressure  gauge,  manometer,  mercury  thermometer, 
speedometer,  voltmeter). 


(3) 


Transducer  - The  desired  measurement  variable  induces  an 
electrical  signal  through  the  transducer  device  which  is 
directly  relatable  to  the  measurement  variable  through 
calibration  (e.g.,  thermocouple,  strain  gage,  photo- 
multiplier, turbine  flowmeter). 


(*) 


Diagnostic  - The  desired  measurement  variable  is  relatable 
to  a direct  measurement  observable  through  an  ab  initio 
model  of  the  physics  (e.g.,  velocity  from  Doppler  frequency 
shift,  temperature  from  quantum  state  population  distribution, 
particle  size  from  scattered  radiation,  molecular  density 
from  Raman  scattered  radiation). 


The  measurements  in  categories  (1),  (2),  and  (3)  are  straightforward,  although 
development  of  sensors,  electrical  signal  conditioning,  and  computer  data 
handling  software  often  presents  formidable  problems.  The  area  of  work 
associated  with  such  problems  is  classed  as  "instrumentation"  development. 
Category  (k) , "diagnostics,"  requires  phenomenological  knowledge  of  the 
fundamental  laws  and  atomic/molecular  properties,  and  thus  often  requires 
a considerable  amount  of  scientific  research.  This  document  is  primarily 
concerned  with  the  assessment  of  the  diagnostics  areas  and  the  required 
fundamental  research  required  to  place  the  various  techniques  on  sound 
footing.  Instrumentation  development  (category  (3))  will  be  discussed 
only  when  the  lack  of  adequate  devices  presents  a barrier  to  the  develop- 
ment of  the  diagnostic  technique. 
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Combustion  flows,  rocket  exhausts  in  particular,  present  a very  hostile 
environment  to  proces  immersed  in  the  flow.  In  addition,  the  presence  of 
the  probe  often  disturbs  the  flow  field  and  the  chemistry  so  that  a repre- 
sentative measurement  cannot  be  made  on  the  gases  sampled.  Consequently, 
immersion  sampling  probes  or  other  sensing  objects  into  the  flow  are 
considered  undesirable,  and  nonlntrusive,  in  situ  techniques  must  be 
employed.  Such  techniques  encompass  emission,  absorption,  and  various 
scattering  mechanisms  of  electromagnetic  radiation. 

AFRPL  interests  in  measurements  are  restricted  to  the  temperature  and 
pressure  ranges  encountered  in  rocket  combustion  chambers,  nozzles,  and 
exhaust  plumes.  In  the  latter  case,  a wide  range  of  temperature  and 
pressure  ranges  are  of  interest,  depending  upon  the  altitude  regime  and 
the  portion  of  the  plume  flow  field  to  be  considered.  The  species  of 
concern  are  the  combustion  products  of  typical  liquid  and  solid  pro- 
pellant rockets.  This  document  is  devoted  to  assessing  the  techniques 
which  have  application  to  rocket  propulsion  rather  than  to  combustion 
diagnos' ics  as  a whole.  Others,  such  as  the  following,  have  addressed 
the  assessment  of  techniques  of  combust  ion  diagnostics  for  other  applications: 


Reference 

Year 

Authors 

Appl ication 

1 

1976 

R.  Goulard,  A.M.  Mel  lor 
and  R.W.  Bilger 

Air  Breathing 

Propulsion  Engines 

2 

1976 

R.  Goulard 

General  Combustion 

Flows 

3 

1977 

S.A.  Self  and 

C.H.  Kruger 

MHD  Combustion 

h 

1977 

A.C.  Ekberth, 

P.A.  Bonczyk  and 

J.F.  Verdieck 

Turbine  Engine  Com- 
bustion and  Flows 

5 

1977 

H.W.  Coleman, 

D.R.  Hardesty,  et.  al. 

Coal  Fired  Gas  Turbines 
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In  addition  to  the  various  reviews,  a number  of  workshops  on  combustion 
diagnostics  have  taken  place  in  the  past  few  years,  at  least  two  of  which 
were  devoted  to  rocket  combustion  diagnostics.0  ^ A valuable  contribution  to 
the  literature  on  combustion  diagnostics  was  recently  made  by  the  publication 
of  a collection  of  papers  given  in  1976  at  the  American  Institute  of  Aero- 
nautics and  Astronautics  (AIAA)  Aerospace  Sciences  Meeting  (Washington,  D.C., 

o 

January  1976).  The  papers  describe  a wide  variety  of  techniques  that  have 
been  applied  to  gas  phase  and  particle  diagnostics  in  combustion  flows. 

Although  the  material  presented  in  References  1 through  8 covers  generally 
the  measurement  problem  areas  applicable  to  rocket  propulsion,  none  of  them 
adequately  brings  into  focus  the  pragmatic  problems  associated  with  "real 
world"  measurements  in  rocket  combustion  flows. 

The  approach  to  be  taken  in  this  document  is  somewhat  different  than  that  of 
other  reviews  and  assessments.  Techniques  which  are  well  established  will 
be  considered  in  the  main  body  of  the  report  only  from  the  standpoint  of 
applicability  and  adequacy  for  rocket-oriented  problems,  and  preferred 
references  will  be  given  to  cite  needed  data  to  support  the  conclusions 
drawn.  The  basic  principles  of  the  established  techniques  will  be  outlined 
briefly  in  the  appendices.  The  plan  of  the  report  is  to  first  establish  the 
needed  measurement  variables,  their  range,  and  the  required  accuracy,  based 
upon  the  end  objective  of  the  measurement  (performance,  plume  signature,  etc.). 
The  available  techniques  will  then  be  listed,  status  of  development  and  ranges 
of  application  will  be  discussed,  and  areas  of  further  development  will  be 
identified  where  applicable.  Of  particular  importance  is  the  identification 
of  voids  in  current  capability  and  introduction  of  possible  new  approaches. 

The  feasibility  of  the  proposed  approaches  will  then  be  studied  in  a pre- 
liminary sense  and  a research  and  development  procedure  suggested. 


2. 


MEASUREMENT  NEEDS  FOR  ROCKET  PROPULSION  FLOWS. 


2. I General  Properties  of  Combustion  Gases. 

The  combustion  process  in  rockets  may  produce  species  in  gas,  liquid,  and 
solid  phases.  The  specific  products  produced  and  their  phase  depends  upon 
the  location  within  the  engine  flow  (combustion  chamber,  nozzle,  or  plume) 
and  the  propellants  being  used.  It  is  not  possible  to  express  many  generali- 
zations without  specifying  the  system  of  interest.  For  the  purposes  of  this 
study,  the  classes  of  propellants  divide  naturally  into  those  producing  only 
gaseous  combustion  products  and  those  producing  an  appreciable  concentration 
of  liquid  or  solid  material.  Cryogenic  propellants  utilizing  liquid  oxygen 
with  liquid  hydrogen  or  hydrazine;  storable  propellants  utilizing  oxidizers 
such  as  nitrous  tetroxide;  fuels  such  as  mono-methyl-hydrazine  or  un symmetrical  - 
dimethyl-hydrazine;  or  any  of  several  other  such  fuels  produce  combustion 
products  which  are  almost  exclusively  gaseous.  Liquid  oxygen  with  kerosene- 
type  fuels  produce  solid  carbon  soot  in  addition  to  the  gaseous  products, 
but  little  attention  is  given  to  hydrocarbon  fuels  in  advanced  propulsion 
systems.  Solid  propellant  systems  may  be  either  highly  metallized,  producing 
large  effluxes  of  liquid  and/or  solid  particulate  matter,  or,  in  the  case  of 
the  smokeless  type  propellants  under  current  development,  may  contain  only  a 
small  concentration  of  particulate  matter.  Advanced  propellants  using  holo- 
genic  oxidizers  and  boron-containing  fuels  produce  gaseous  and  liquid/solid 
products,  as  yet  not  well  defined,  and  little  attention  will  be  given  to  this 
class  of  propellants  in  this  study. 

Appendix  A contains  tables  of  properties  and  species  concentrations  for 
several  propellant  combinations  based  upon  one-dimensional  equilibrium  cal- 
culations. The  levels  of  species  concentrations  and  temperatures  listed  there 
serve  as  references  for  this  study.  The  major  species  produced  in  the  combus- 
tion processes  are: 
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Propel lants 

Vj 

vc*h, 

WVyHz 

Sol  ids 


Major  Gaseous  Species 

h2o,  h2 
h2o,  co2,  CO 

N2'H2°*  C02’  H2*  C0 
N2,H20,  Hc£,  C0,H2 


Major  Solid  Species 


A£2°3 


The  gaseous  species  listed  generally  make  up  about  99  percent  of  the  total 
gas,  with  the  remainder  composed  of  poorly  known  concentrations  of  minor 
stable  species,  such  as  unburned  fuel  fragments,  oxides  of  nitrogen,  alkali 
metals  (Na  or  K) , and  free  radicals  (0,  OH,  CH,  etc.).  Thus,  gas  phase 
diagnostics  is  directed  toward  one  or  more  of  the  major  molecules,  unless 
a specific  need  exists  to  determine  concentration  of  one  of  the  minor  species. 

Solid  materials  produced  consist  of  either  carbon  soot  (C  ) for  the  0,/C  H 

n z x y 

liquid  propellant  and  some  advanced  propellants,  such  as  C£F^/MMH,  or  metal 
oxides  (mostly  AJljO^)  for  the  solid  propellants.  Thus  diagnostics  of  partic- 
ulates in  rocket  combustion  should  be  concerned  mostly  with  carbon  soot  or 
AHjOj  particles.  Solid  propellant  rockets  which  do  not  use  metals  imbedded 
within  the  propellant  grains  do  not  produce  an  appreciable  number  of  partic- 
ulates in  the  combustion  process,  but  condensation  In  the  far  field  of  the  low 
altitude  plume  often  occurs.  The  condensate  (secondary  smoke)  is  thought  to 

be  composed  of  H„0/HC£  droplets  centered  on  small  nucleation  sites  with  the 

q 

possibility  that  HC£  plays  a role  in  the  condensation  process.  Such 
condensates  are  also  subject  to  diagnostics. 

In  Figure  1,  a schematic  of  a combustion  flow  for  a liquid  propellant  rocket 
is  shown  for  a low  altitude  afterburning  plume  (Figure  la)  and  a vacuum  plume 
(Figure  lb).  In  Figure  2,  a similar  schematic  is  given  for  a solid  propellant 
rocket  flow  at  low  altitude  (air-launched).  A temperature  and  pressure  history 
along  the  flow  paths  from  the  combustion  chamber  to  the  far  field  plume  is  given 
in  each  case.  These  schematics  will  be  used  as  maps  to  indicate  the  applica- 
bility of  different  diagnostic  methods.  Those  methods  which  have  been  used 
are  shown  at  the  appropriate  location  In  the  flow  schematic.  Figures  1 and  4 


A 


-II- 


PRIMARY  SMOKE  REGION  NO  VISIBILITY  SECONDARY  SMOKE  REGION 


will  be  discussed  in  more  detail  in  the  next  section.  The  intent  at  this 
point  is  to  define  the  ranges  of  variables  to  be  dealt  with  as  indicated 
on  the  figures. 


j i 
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2.2  General  Needs  for  Measurements  in  Combustion  Gases. 


2.2.1  Performance  of  Liquid  Propellant  Rockets. 


In  the  development  of  rocket  propulsion  systems  to  meet  the  requirements  of 
users  like  the  Space  and  Missile  Systems  Organization  and  Aeronautical  Systems 
Division,  the  major  requirement  is  performance,  usually  measured  in  terms  of 
"specific  impulse"  (I  ),  which  is  defined  as 


I 


sp 


Thrust  (F) 

Weight  flow  of  propellants  (W) 


sec 


(1) 


The  thrust  and  weight  flow  of  propellants  in  a liquid  propellant  rocket  are 

directly  measurable  so  that  I is  instantaneously  determined  throughout  a 

firing.  The  energy  release  by  the  chemical  reaction  of  fuel  and  oxidizer  is 

known  from  chemical  heat  of  formation  data  so  that  the  ideal  I,  can  be 

sp 

predicted.  A measure  of  the  performance  is  then  available  through  comparison 
of  the  measured  to  the  ideal  I . In  the  actual  case,  the  chemical  reactions 
rarely  go  to  completion,  and  this,  coupled  with  real  gas  dynamic  behavior, 
prevents  the  achievement  of  perfect  performance.  The  measurement  of  F and 
W in  a gross  sense  does  not  indicate  the  reason  for  a less  than  ideal  per- 
formance; thus,  because  these  have  been  traditionally  the  only  measurements 
made,  development  of  rocket  engines  has  evolved  as  an  empirical  cut-and-try 
procedure. 


With  the  advent  of  the  modern,  high-speed  digital  computer  and  improved 
understanding  of  turbulent,  chemically  reacting  flow  and  droplet  burning 
phenomena,  attempts  have  recently  been  made  at  putting  the  development  of 
liquid  propellant  rockets  on  a sound  scientific  basis  through  analytical 
modeling  (Reference  10,  for  example).  However,  such  complex  models  involve 
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many  assumptions,  and  gross  measurements  of  F and  W are  seldom  sufficient 
to  determine  just  where  a model  is  deficient.  For  example,  thrust  F is 
determined  by 


where 


F - 5 ' / * l - F0)  CO,  e 2„rw  dz  (2) 

g is  the  acceleration  of  gravity; 

p is  the  gas  density  at  some  position  x , y in  the  plume  of  the 
G c c 

nozzle  exit; 

V is  the  gas  velocity  at  some  position  x , y in  the  plane  of  the 
e e e 

nozzle  exit; 

dA  is  a differential  area  at  a position  x , y in  the  plane  of 
0 “ “ 

the  nozzle  exit; 

is  the  static  pressure  at  the  wall  at  an  axial  position  z from 
the  nozzle  throat; 

P is  the  static  pressure  of  the  ambient  surroundings; 
o 

0 is  the  angle  which  the  nozzle  makes  with  the  z axis  of  the  nozzle 

r is  the  radius  of  the  nozzle  at  the  axial  position  Z; 
w 

L is  the  axial  length  of  the  nozzle  from  the  throat  to  the  exit 
plane. 


The  gas  density  and  velocity  at  any  point  in  the  exit  plane  are  determined 
by  the  oxidizer-to-fuel  ratio  (0/F)  at  that  point  and  the  degree  to  which 
combustion  has  gone  to  completion,  as  well  as  the  gas  dynamic  conditions 


imposed  by  the  nozzle  expansion.  Thus,  an  accurate  computer  program  for 
predicting  the  I must  start  at  the  propellant  injector  face  and  follow 
the  mixing  and  combustion  through  the  combustion  chamber  to  the  nozzle 
exit.  Then  it  must  compute  the  integrals  of  Equation  (2)  point  by  point. 

A measurement  to  test  the  adequacy  of  the  program  must  include  the  local 
exit  plane  density  and  velocity  or  sufficient  other  variables  to  separately 
calculate  them  so  that  point  by  point  comparison  between  model  prediction 
and  measurements  can  be  made.  A sufficient  set  of  variables  to  compute 
Equation  (I)  are  the  temperature.  Tfi,  the  0/F,  the  total  pressure  profiles, 
and  a mapping  of  P (z).  The  pressures  are  directly  measurable,  but  T and 
0/F  are  not;  this  is  the  beginning  of  the  diagnostics  problem.  Measurement 
of  the  temperature  and  the  major  species  concentrations  (H^O,  CO^,  CO,  and 
perhaps  or  Hg)  is  sufficient  to  determine  a local  value  of  the  0/F.  If 
certain  assumptions  are  made  about  the  chemistry,  then  measurement  of  the 
local  density  of  only  one  or  two  of  these  species  may  be  sufficient. 

2.2.2  Performance  of  Solid  Propellant  Rockets. 

In  the  testing  of  solid  propellant  rockets,  the  weight  flow  rate  of  propel- 
lants,  W,  is  not  directly  measurable  but  depends  upon  the  propellant  surface 
burning  rate.  Since  the  burning  rate  is  also  not  measurable,  the  ls  is 
usually  determined  by 


I 

sp 


Fdt 


(3) 


i 


where  l(  is  the  total  impulse; 

W is  the  total  weight  of  propellant; 
tQ  is  the  total  burn  time. 
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A model  to  compute  the  thrust  F for  solid  propellant  rockets  follows 
Equation  (2)  except  that  now  both  gas  and  solid  phases  must  be  taken  into 
account  In  computing  P , Vg,  and  Py.  Consequently,  the  temperature  of  both 
phases  must  be  computed,  and  the  composition  of  both  gas  and  solids  must  be 
determined.  The  interaction  of  gas  and  particle  flow  is  a major  part  of 
the  analytical  modeling  for  solid  propellant  rocket  motors,  and  a critical 
variable  is  the  status  of  the  particles  in  the  flow  (size,  spatial  distri- 
bution, temperature,  and  constituency  of  the  particles).  By  constituency 
is  meant  the  chemical  composition,  the  porosity,  the  geometrical  shape,  and 
the  phase  (liquid,  solid,  crystalline,  etc.)  of  the  particles.  The  only 
information  currently  available  on  particle  properties  is  that  obtained 
through  analysis  of  samples  obtained  during  or  after  a motor  firing.  Such 
particle  analyses  have  shown  that  carbon  soot  with  a mean  particle  diameter 
well  below  1 ym  is  formed  in  the  combustion  of  hydrocarbon  fuels  (see,  for 
example,  Reference  11). 

For  solid  propellant  motors  to  which  aluminum  has  been  added,  the  particles 

collected  are  AJt-0,,  with  a wide  range  of  size  distributions,  the  size 

£ } 12 
depending  to  a large  extent  on  how  the  samples  were  obtained.  In  some 

cases,  the  AJl^O^  Particles  are  porous,  or  even  hollow,  with  evidence  that 

they  were  in  a liquid  state  within  the  nozzle. 

The  temperature  history  of  the  particles  depends  upon  their  history  in  the 
flow  and  on  their  radiative  properties.  The  temperature  of  the  particles 
may  not  be  the  same  as  the  gas  temperature  because  of  the  gas/particle  heat 
and  mass  and  state  changes. 

The  measurements  required  to  fully  evaluate  analytical  models  for  the 
performance  of  solid  propellant  rockets  are  summarized  as  follows: 

Gas  Phase:  Local  gas  temperature 

Local  composition  of  major  species 
Velocity 
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Condensed  Phase: 


Geometrical  size  distribution 


Composition 
Hass  distribution 
Velocity  distribution 
State  (liquid  or  solid) 

Temperature 

Hardly  any  of  these  variables  is  directly  measurable  by  probes  because  of 
disturbances  induced  at  the  probe  entrance.  Consequently,  in  situ  diagnostic 
methods  are  required.  The  diagnostic  problem  for  the  solid  propellant  rocket 
is  clearly  much  more  formidable  than  for  the  liquid  propellant  rocket. 

2.2.3  Rocket  Plume  Signatures. 

The  AFRPL  is  also  responsible  for  providing  information  to  users  on  the 
electromagnetic  observables  and  interferences  produced  by  rocket  exhaust 
plume  products.  The  observables  of  interest  may  lie  in  any  wavelength 
range  from  vacuum  ultraviolet  to  radio  frequency.  The  applications  include 
surveillance  of  enemy  missile  launches,  missile  warning  receivers  for  air- 
craft, radar  cross  sections,  transmission  or  communications  signals  through 
plumes,  contamination  of  spacecraft  surfaces  by  plume  materials,  and' various 
countermeasures  to  disguise  or  hide  missile  launches  and  maneuvers.  The 
extent  of  the  flight  regimes  and  scenarios  is  much  too  broad  to  rely  on 
flight  test  or  ground  test  programs  to  supply  complete  information  over  all 
flight  regimes.  Therefore,  reliance  is  made  on  analytical  gas  dynamic, 
chemistry,  and  electromagnetic  radiation  models  to  provide  the  overall  pre- 
dictive capability.  Two  measurement  problems  arise:  (a)  The  models  require 
input  data  (initial  conditions),  and  (b)  they  must  be  validated  through 
appropriate  testing  to  insure  that  the  most  Important  physical  processes 
have  been  modeled. 


-19- 


In  order  to  illustrate  the  measurements  required,  consider  Figure  3 in  which 
a plume  is  shown  which  is  emitting  radiation  and  may  also  scatter  radiation 
from  other  sources  to  an  observer  (detector).  Depending  upon  the  orienta- 
tion of  the  detector  with  respect  to  the  plume,  the  wavelength  range,  and 
the  orientation  of  the  external  radiation  source.  Figure  3 represents  a wide 
range  of  scenarios.  A generalized  expression  for  the  radiation  received  at 
the  detector  may  be  written  as  follows: 


'L  ■ 'm  Ui  hl  h T <T-V»>  d'dC!d> 

9t(T,P.,X) 

+ ;ax  ;AZ  ;Ae  'o  (9’x)  sisx  (0’pi'x) dzdfldXde  (M 


where  is  the  detected  radiation  within  some  specified  wavelength 

range  AX  * X^  - X^; 

E(T,P.,X)  is  the  emission  coefficient  of  the  radiation  at  some 
position  with  respect  to  the  detector  produced  by  species  i at 
temperature  T and  partial  pressure  P.  at  a wavelength  X; 

t(T,P.,X)  is  transmittance  of  the  radiation  at  a wavelength  X along 
the  path  z to  the  detector; 

l.(B>X)  is  the  irradiance  from  an  external  source  falling  on  the 
plume  at  an  angle  6 with  respect  to  z (such  as  the  sun  or  a radar 
beam)  v,.'  may  be  scattered  into  the  detector  field  of  view  or 
be  attenuated  by  the  plume; 

S.(e,P.,X)  is  the  scattering  function  of  the  plume  at  the  angle  P 
by  species  (particles)  i at  the  wavelength  X. 

I 
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The  integrals  are  carried  out  over  the  wavelength  of  interest  (the  bandpass 
of  the  detector),  the  possible  directions  of  the  external  radiation,  the 
solid  angle  subtended  by  the  detector,  and  the  distance  along  the  vector  z. 

The  reason  for  writing  the  transfer  equation  with  this  much  generality  was 
to  show  the  commonality  of  the  requirements  for  knowledge  of  the  temperature 
and  species  densities  throughout  the  plume  for  all  predictions  of  electro- 
magnetic signature  properties.  Thus,  a prerequisite  to  a radiation  model  is 
a gas  dynamic-chemistry  model  capable  of  accurate  prediction  of  local  tempera- 
ture, local  gas  species  density,  and  particle  size  distribution  and  thermo- 
physical state. 

The  plume  signature  and  plume-electromagnetic  wave  interaction  modeling  area 
actually  encompasses  a wide  spectrum  of  separately  definable  problems. 

Some  of  the  major  modeling  requirements  are  listed  as  follows: 

(1)  Prediction  of  (ly  to  25y)  radiant  intensity  from  liquid  and 
solid  propellant  strategic  and  tactical  missile  plumes. 

a.  Low  altitude  (0  to  40  km)  afterburning  plumes. 

b.  High  altitude  (~50  to  70  km)  vacuum  plumes. 

(2)  Prediction  of  surface  contaminant  species  from  space  orbit- 
control  thrustors. 

(3)  Prediction  of  radar  cross  section  of  missile-plume  combination. 

(4)  Prediction  of  microwave  signal  attenuation  and  scattering  upon 
passage  through  tactical  missile  plumes. 

All  of  these  models  require  as  input  data  the  exit  plane  profiles  of  static 
temperature  and  partial  pressures  of  major  species,  the  particulate  properties 
(if  present),  and  a total  pressure  exit  plane  profile.  Once  validated  through 
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experiments  employing  proper  diagnostic  measurements,  the  performance  models 
discussed  in  2.2.1  and  2.2.2  can  provide  the  exit  plane  inputs  to  the  plume 
models.  Evaluation  of  the  plume  signature  prediction  models  can  be  accom- 
plished in  an  "at  large"  sense  by  measurement  of  the  predicted  variable  in 
either  flight  test  or  captive  simulation  tests.  However,  field  tests  are 
seldom  satisfactory,  and  the  deficient  elements  of  the  model  cannot  be 
detected.  Therefore,  controlled  simulation  testing  is  necessary  in  which 
measurements  of  critical  variables  (temperature,  major  species  concentration, 
electron  density,  etc.)  at  selected  positions  in  the  plume  are  made. 

Other  plume-produced  phenomena  of  interest  include  the  following: 

(5)  The  interaction  of  plume  species  with  atmospheric  species  to 
produce  radiation  in  the  far  field  of  high  altitude  plumes. 

(6)  The  condensation  of  particulates  in  the  far  field  of  smokeless 
propellant  tactical  missile  plumes  to  form  condensation  trails. 

These  phenomena  require  a different  kind  of  diagnostics  than  the  near  field 
plume  phenomena.  The  high  altitude  interaction  problem  requires  measurements 
of  minor  or  trace  species  concentrations  at  low  density.  The  low  altitude 
condensation  problem  requires  particle  property  diagnostics  of  the  condensates. 
Neither  of  these  phenomena  can  be  studied  adequately  by  flight  testing  or  by 
rocket  firings  in  captive  facilities,  but  require  simulation  of  the  phenomena 
in  some  indirect  way,  making  good  diagnostics  that  much  more  important. 

2.2.  b Fundamental  Mechanism  Studies. 

Many  fundamental  studies  are  required  in  the  development  of  new  concepts  in 
rocket  propulsion  and  in  understanding  of  plume  phenomena.  The  laboratory 
experiments  conducted  during  these  studies  also  require  diagnostics.  Some 
fundamental  studies  which  currently  can  be  identified  are  listed  below, 
along  with  the  measurement  desired: 
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(1)  Solid  propellant  combustion. 


Objective:  To  determine  burning  rate  and  surface  reaction 
phenomena  for  a variety  of  propellants. 

Apparatus:  Combustion  bomb. 

Measurement:  Geometrical  history  of  reacting  surfaces  and 
escaping  particles  (with  definition  of  gaseous  species 
properties  if  possible). 

(2)  Metal  combustion  phenomena. 


Objective:  To  determine  the  mechanism  of  metal  oxide  forma- 
tion, nucleation,  and  particle  growth. 

Apparatus:  Low  density  reaction  furnace. 

Measurement:  Gaseous  specie  identification,  concentration, 
and  temperature  and  particle  growth  history  and  constituency. 

(3)  UV  radiative  mechanisms. 

Objective:  To  identify  chemiluminescent  excitation  mechanisms 
for  UV  plume  radiators  (such  as  CO,  NO,  OH,  CN,  CH,  and 
continuum  mechanisms)  and  determine  rate  constants. 

Apparatus:  Gas  reactor  configuration,  shock  tubes,  etc. 

Measurement:  Species  concentrations  and  temperature  of 
reaction  participants,  and  UV  radiation. 
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(4)  Particle  radiative  properties. 


Objective:  To  determine  spectral  emissivity  and  scattering 
cross  sections  of  particles  as  encountered  in  combustion 
gases. 

Apparatus:  Two-phase,  heated  flow  facility. 

Measurements:  Radiation  and  scattering  function  and  temper- 
ature. 

(5)  Spectral  properties  of  gas  species. 

Objective:  To  determine  the  radiative  constants  (molecular 
properties  or  band  model  properties)  for  molecules  of 
interest  in  plumes. 

Apparatus:  Burners  and  heated  flow  facilities. 

Measurements:  Temperature,  species  density,  and  spectral 
radiation. 

These  are  but  a few  of  the  many  fundamental  programs  being  conducted  by  the 
AFRPL  which  require  diagnostic  type  measurements  for  proper  treatment  of 
the  experiments. 

2. 3 Summary  of  Diagnostic  Needs. 

The  measurement  needs  discussed  in  2.2.1,  2.2.2,  and  2.2.3  are  summarized 
in  Tables  la  and  lb  for  liquid  propellant  and  solid  propellant  rockets, 
respectively.  The  variable  to  be  measured,  the  region  in  the  plume  where 
the  measurement  is  desired,  the  applications  involved,  the  ranges  of  the 
variables  of  interest,  and  the  accuracy  required  to  accomplish  the  objective 
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TABLE  lb.  SUMMARY  OF  MEASUREMENT  NEEDS  FOR  ROCKET  PROPULSION  APPLICATION  - SOLID  PROPELLANTS  (1  of  2). 
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of  the  application  are  given  in  the  tables.  The  information  contained  in 
Tables  la  and  lb  will  be  used  to  assess  the  capability  of  existing  diagnostic 
techniques  and  to  locate  inadequacies  so  that  new  work  can  be  planned. 


! 
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The  emission  or  absorption  of  radiation  from  atoms  or  molecules  in  the  gas 
phase  has  the  property  that  the  absolute  value  is  related  directly  to  the 
number  density  of  the  specie,  and  the  spectral  distribution  is  a function 
of  the  temperature  through  some  formulation  of  the  quantum  state  equilibrium 
distribution.  Several  procedures  have  been  used  to  obtain  a "measure"  of 
the  absolute  distribution  function  of  a specie  in  a gas  flow,  including 


Passive  emitted  radiation 


Absorbed  radiation  from  a source 
Fluorescent  radiation  induced  by  a source 
Fluorescent  radiation  induced  by  electron  beams 
Scattering  from  an  intense  radiation  beam  (laser)  at  frequency  vq 
Rayleigh  - nonshifted,  vq 
Raman  - shifted,  v + vn,  vD  * Raman  shift 

O ™ r>  K 

Nonlinear  Raman 


From  these  phenomena  have  come  many  specialized  techniques  for  gas  analysis 
and  measurement  of  properties  of  ionized  gases  (plasmas),  flames,  and  combus- 
tion gas  flows  from  turbojet  and  rocket  propulsion  systems.  The  only  tech- 
niques to  be  discussed  here  are  those  which  have  been  applied  to  combustion 
flows  and  those  for  which  there  appears  to  be  no  barriers  to  use  in  combustion 
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flows  since  they  have  been  used  successfully  in  similar  flows.  The  tech- 
niques are  listed  in  Table  2 along  with  the  measurements  possible  with  them 
and  the  ranges  of  application.  The  physical  principle,  instrument  layout, 
ranges  of  applicability,  and,  where  possible,  the  estimated  typical  uncer- 
tainty of  each  of  the  techniques  listed  in  Table  2 are  discussed  in 
Appendix  B. 

The  measurement  of  velocity  in  gas  flows  by  an  in  situ  technique  depends 
upon  use  of  the  Doppler  frequency  shift  of  an  emission  spectral  line  or  on 
the  Doppler  shift  of  scattered  radiation  from  a laser  beam.  In  practical 
combustion  flows,  the  Doppler  shifts  of  emission  lines  are  seldom  large 
enough  relative  to  the  line  width  due  to  random  motion  and  collisional 
broadening  to  be  measurable.  The  technique  that  has  proved  most  successful 
for  velocity  measurement  has  been  the  Laser  Doppler  Velocimeter  (LDV) 
technique  in  which  the  frequency  shift  of  the  scattering  signal  from 
particles  imbedded  in  the  gas  flow  is  measured.  (This  technique  is 
discussed  at  length  in  the  survey  References  1 through  5 and  by  several 
papers  in  Reference  8).  The  discussion  of  this  technique  for  rocket  pro- 
pulsion applications  properly  belongs  in  the  next  section  on  gas/particle 
flows.  Quite  obviously,  the  chief  question  with  regard  to  measurement  of 
gas  flow  velocities  has  to  do  with  whether  particles  are  moving  at  the  same 
velocity  as  the  gas  flow.  Generally,  particles  having  diameters  less  than 
about  0.5  ym  will  follow  the  gas  flows  of  rockets  at  pressures  greater  than 
about  1 atm.  The  question,  however,  is  a much  more  complex  one  which 
requires  careful  analysis  for  each  application  (see  References  31*  and  35, 
for  example). 

The  principal  application  of  diagnostic  methods  to  liquid  rocket  combustion 
flows  has  been  the  measurement  of  nozzle  exit  plane  temperature  and  species 
concentration  profiles  in  simulated  vacuum  or  flight  environments.  Thus, 
the  discussion  to  follow  will  center  mostly  on  these  applications. 
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The  emission  and  absorption  of  infrared  radiation  by  the  principal  products 

of  combustion  of  conventional  liquid  propellants,  t^O  and  CC^,  provides  a 

means  of  measurement  of  temperature  and  the  partial  pressure  of  H^O  and  COj 

in  the  exhaust  of  rocket  engines.  The  IR  E/A  technique  is  based  upon  the 

1 3 

band  model  radiative  transfer  concept  and  has  been  used  on  a number  of 

applications  to  rocket  test  ng.  (See,  for  example.  Reference  14.)  Two  of 

the  most  recent  applications  are  given  in  References  15  and  16,  where  it  is 

shown  that  small  changes  in  propellant  injector  configurations  are  detectable 

15 

using  the  technique.  In  one  application,  the  results  of  applying  the 
technique  have  been  used  to  make  a preliminary  assessment  of  liquid  engine 
performance  computer  codes.  Current  effort  on  the  technique  is  being 
devoted  to  improvements  in  instrumentation  and  data  treatment  with  the  goal 
being  to  provide  higher  confidence  in  the  data  through  better  uncertainty 
analysis. Current  testing  efforts  are  devoted  to  extensions  of  the  per- 
formance code  evaluation  as  described  in  Reference  16  in  order  to  provide 

1 8 

validation  for  the  latest  performance  codes. 


Resonance  Absorption 


sndix  B-2) . 


The  ultraviolet  (UV)  resonance  absorption  method  is  applicable  to  minor 
diatomic  species  and  free  radials  such  as  NO,  OH,  NH,  CN,  and  CH.  Applica- 
tion to  rocket  combustion  flows  has  not  been  made,  but  results  of  several 

I8-20 

measurements  in  turbojet  engine  exhausts  have  been  published.  Labora- 

tory equipment  has  so  far  been  employed  for  this  technique.  However,  design 
of  a dedicated  system  would  not  be  difficult.  Estimates  of  uncertainty 
(_  + 15  percent)  are  probably  conservative,  and  a thorough  analysis  has  not 
been  done.  No  need  for  application  to  rocket  exhaust  flows  is  anticipated 
at  this  time. 


3. i.1.3  Rayleigh/ Raman  Scattering 


indices  B-3  and  B-4) 


22 

Rayleigh  scattering  as  a measure  of  gas  density  and  Raman  scattering  as 

24 

a measurement  technique  for  temperature  and  species  concentration  have 
received  much  attention  in  the  past  few  years. 
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Rayleigh  and  Raman  scattering  data  from  a rocket  exhaust  were  first  obtained 

during  an  AFRPL  sponsored  test  program  in  I97W5.  The  program  served  as 

a feasibility  demonstration,  and  temperature  and  species  concentration  data 

were  successfully  obtained.  Many  problems  (background  interference, 

fluorescence,  plume  fluctuations,  etc.)  were  encountered  in  the  test,  and 

the  experience  gained  was  used  in  designing  apparatus  for  use  of  Raman 

25 

spectroscopy  in  a subsequent  test  program.  In  the  latter  test,  the 
objective  was  to  measure  instantaneous  temperature  and  number  density  of 

in  a turbulent  mixing  zone  between  external  air  flow  (Mach  number  2 and  3) 
and  a rocket  plume.  A laser  pulse  width  of  about  20  nanoseconds  was  used  at 
a repetition  rate  of  12  pulses  per  minute  using  a Q-switched  ruby  laser. 

Oata  uncertainty  was  about  +_  8 percent;  fluctuations  in  temperature  up  to 
+_  20  percent  were  measured.  The  limitations  imposed  in  both  tests  were 
mainly  economic  in  that  better  data  requires  either  much  more  instrumentation 
or  much  longer  test  times.  It  is  concluded  that  the  Raman  scattering  tech- 
nique can  be  used  to  gather  most  data  required  in  gaseous  combustion  flows, 
provided  funds  are  available.  Emphasis  is  currently  being  placed  on  improve- 
ments in  equipment,  such  as  the  use  of  image  intensifiers  to  obtain  profile 
data  and  Fabry-Perot  etalons  for  higher  throughput  than  can  be  obtained  using 
spectrometers. 

Nonlinear  Raman  excitation  techniques,  such  as  Coherent  Anti-Stokes  Raman 

2 6 

Spectroscopy  (CARS),  are  still  in  the  development  stage  and  were  not 
considered  in  this  assessment. 

3 . 1 . 1 . 4 High  Resolution  Spectroscopy  (Appendix  B-5). 

Molecules  which  have  a well  defined  vibrat ional -rotat ional  line  structure, 
such  as  CO  or  HF,  are  approachable  either  in  emission  or  absorption  by 
analysis  of  the  resolved  line  structure.  HF  has  been  used  successfully  as 
a tracer  gas  for  mixing  studies  by  absorption  spectroscopy.  The  applica- 
tion to  rocket  gases  using  conventional  high  resolution  spectroscopy  is 
discussed  in  Reference  28,  where  a high  resolution  CO  spectrum  from  a rocket 


exhaust  is  shown,  and  the  criteria  for  diagnostics  of  CO  are  defined.  The 

application  of  tunable  lasers  to  make  high  resolution  absorption  measurements 

29 

in  flames  is  also  a very  attractive  technique.  The  difficulty  with  the 
high  resolution  method  is  increased  over  the  IR  band  model  method  by  the 
necessity  to  solve  the  inversion  problem  for  the  radiative  transfer  on  a 
line  by  line  basis.  However,  the  end  result  using  resolved  lines  will  have 
a higher  confidence  level  than  for  the  band  model  because  the  physics  is 
much  better  defined.  The  conventional  spectroscopy  is  straightforward,  and 
no  substantial  barriers  to  application  are  anticipated  for  molecules  having 
a resolvable  line  structure. 


3. 1.1.5  Electron  Beam  Fluorescence  (Appendix  B-6) . 

Electron-beam  excitation  of  fluorescent  radiation  can  be  used  to  provide 

a measure  of  the  distribution  function  for  several  molecules  and  atoms  of 

interest  in  far  field  rocket  plumes  under  high  altitude  simulation  conditions. 

The  technique  is  useful  in  a direct  sense  for  gas  number  densities  of  about 

lo'^  cm  ^ (the  signal  sensitivity  limit)  to  about  10*^  cm  ^ (upper  limit 

governed  by  quenching  collisions).  A typical  example  is  given  by  the 

application  to  chemical  laser  type  flows^  in  which  species  concentration 

measurements  in  the  mixing  region  between  adjacent  jets  were  made.  The 

latter  measurements  also  demonstrated  the  problems  which  can  arise  when 

proper  attention  is  not  given  to  the  complex  energy  transfer  and  quenching 

reactions  that  may  occur.  The  technique  has  been  used  with  considerable 

success  in  low  density  expansions  of  gases  in  which  the  condensation 

31 

properties  were  being  studied. 

3 . 1 . 1 . 6 Resonance  Fluorescence  (Appendix  B-7) . 


In  the  resonance  absorption  technique,  the  transmittance  of  source  radiation 
tuned  to  the  resonance  of  the  absorbing  media  is  measured.  The  source  may 
be  a gas  discharge  tube  containing  the  same  gas  for  which  absorption  measure 
ments  are  desired,  or  it  may  be  a tunable  laser.  The  excited  state  which 
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results  from  absorption  of  a photon  will  either  relax  by  emission  of  a 

photon  (fluorescence)  or  by  collision  with  other  atoms  or  molecules 

(quenching).  If  the  quenching  rate  is  low,  the  emitted  radiation  may  be 

analyzed  in  much  the  same  way  as  E-beam  exalted  fluorescence.  However, 

the  cross  section  for  E-beam  excitation  is  .arger  than  for  fluorescence 

so  that  no  advantage  is  generally  gained  by  using  fluorescence  in  the 

density  range  below  which  quenching  is  a problem.  However,  it  was  recently 

32 

pointed  out  by  Daily  that  in  the  limit  of  a large  amount  of  quenching, 

the  need  for  detailed  knowledge  of  quenching  cross  section  disappears,  and 

the  technique  may  be  used  for  several  species  including  the  alkali  metals. 

The  method,  which  employs  a tunable  laser  in  the  UV  or  Visible,  is  then 

known  as  "saturated  fluorescence."  A successful  extension  of  the  method 

as  outlined  by  Daily  was  made  to  C£  in  flames  in  which  saturation  could  not 

be  obtained,  by  noting  the  linear  dependence  of  the  fluorescence  on  the 
33 

source  power.  The  technique  appears  to  have  potential  for  use  in  rocket 
combustion  flows  containing  a considerable  amount  of  free  carbon. 

3.1.2  Assessment  of  Gas  Phase  Only  Diagnostics. 

From  consideration  cf  the  "needs"  for  liquid  propellant  rocket  combustion 
flow  measurements  (Section  2.1)  and  the  techniques  that  are  available  to 
make  such  measurements  (Table  2),  it  is  concluded  that: 

(1)  Existing  techniques  are  adequate,  in  principle,  to  meet 
currently  anticipated  measurement  needs. 

(2)  Refinements  in  hardware,  data  handling  software,  and  uncer- 
tainty analysis  will  continue  to  require  attention. 

In  other  words,  it  is  concluded  that  basic  research  on  new  techniques  for 
gas  (only)  flows  is  not  required,  while  continued  development  is  required 
to  make  measurements  more  economical  either  from  the  hardware  or  test 
complexity  standpoint,  and  to  improve  sensitivity,  accuracy  and  confidence 
level  of  the  measurements. 
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A discussion  of  the  cost  of  making  a measurement  seems  to  be  in  order.  A 
rough  attempt  to  classify  the  techniques  according  to  the  cost  of  making 
one  radial  profile  measurement  of  the  density  and  temperature  of  a specie 
at  nozzle  exit  conditions  is  given  in  Table  3.  Note,  however,  that  there 
is  little  overlap  in  species  which  can  be  measured  by  more  than  one  tech- 
nique; i.e.,  H20  can  only  be  measured  by  IR  E/A,  C02  can  be  measured  by 
IR  E/A  or  Raman,  CO  can  be  measured  by  either  Raman  or  high  resolution, 

N2  can  only  be  measured  by  Raman,  etc.  Thus,  total  reliance  on  only  one 
technique  for  an  application  is  not  always  possible. 

3. 2 Techniques  for  Gas/Particle  Flows. 

3.2.1  Assessment  of  Gas/Particle  Flow  Diagnostics. 

A survey  of  diagnostic  applications  to  the  gas/particle  flows  produced  by 
solid  propellant  combustion  revealed  no  successful  application  of  either 
particle  or  gas  phase  techniques.  Much  research  on  particle  diagnostics 
by  optical  methods  in  simpler  flow  regimes  (flames,  wind  tunnels,  etc.)  is 
in  progress,  but  the  only  attempts  at  application  to  solid  rocket  plumes 

have  been  unsuccessful.  Such  attempts  include  Laser  Doppler  Velocimetry 

* 36  37 

for  velocity^  and  holography  for  particle  size.  Attempts  at  using  the 

IR  E/A  technique  for  gas  phase  diagnostics  at  the  exit  of  a solid  propellant 

38 

motor  nozzle  likewise  met  with  failure,  and  interpretation  of  spectral 

00 

emission  data  in  terms  of  particle  temperature  has  also  been  unsuccessful. 
Therefore,  it  can  be  concluded  from  the  outset  that  existing  techniques  are 
inadequate  to  meet  the  needs  for  measurements  in  the  gas/particle  flows  of 
sol  id  rocket  motors. 

Interest  in  gas/particle  flows  does  not  arise  solely  from  solid  propellant 
rockets  but  also  from  problems  encountered  with  such  processes  as  coal  fired 
turbine  and  magnetohydrodynamic  power  generation.^  ^ However,  no  actual 
diagnostic  measurements  have  been  reported — only  studies  to  identify 
approaches.  It  is  instructive  to  list  the  approaches  utilized  to  date  for 
in  situ  particle  diagnostic  measurements  as  given  in  Table  k.  The  most 
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TABLE  3.  COST  ANALYSIS  FOR  GAS  PHASE  COMBUSTION  DIAGNOSTIC  TECHNIQUES. 


«/) 

Si 


o — 


UJ  o 


X o 

3 o 


cc  * 
< <0-1 

o h- 

OC  CO 

< o 

X o 


LA 


3 

X 


•fc 

O 


* 

o 


* 

«fc 

* 

o 


3 

O 

X 

(9 


V) 

“D 

C • 
O Li 

o </> 

4)  O 

«0  O 


I 

IA 


* 

* 

«*c 

O 

O 


3 

X 


O 4) 

— «/3 

C X 
<D 

X 4) 


03 

C 


10 
oi  cr» 


4)  g 


•fc 

«fc 

<M 


* 

* 

SO 


O 

sO 


IA 


i 

o 

IA 


O *-» 


03 

c 


(/> 

<V 

O 


to 

tO  L. 

o>  O 

C • ■ 

— u 

l.  4) 


<o 

4> 

>- 


to 

4) 

>- 


to 

4) 

>- 


to 

03 

c 


M-  4) 
“O 

*o 

c U- 
o o 
u 

4)  U 
to  4) 
X 
O E 

sO  3 
C 
4) 

X 4) 


z 

© 

u 

<M 

o 

U 

•k 

•k 

X 

Ok 

o 

u 

u 

IV 

o 

X 

N 

z 

«/> 

O 

«k 

«k 

• 

Ui 

u 

i 

«k 

<s 

u 

u. 

* 

o 

X 

X 

u 

•k 

Ui 

o 

«k 

•k 

•k 

• 

•k 

•l 

N 

o 

X 

U- 

o 

40 

X 

o 

X 

X 

O 

4) 

X 

E 

3 

C 


4) 

O 

3 

"O 

4) 

u. 


tO  Li 
4)  10 
3 O 
CT  U 


l_  4-1 

O <o 
♦-»  o 
u u 

4) 

4->  4) 
4)  L. 
*0  «0 
5 


O 

to 

4) 


03 


4-*  O 


u 

4) 

CO 

<V 


4> 

X 

IV 

C 


4) 

O >* 
D — 
"O  4> 
4>  «M 
L-  (0 
E 

c — 
(V  X 

o o 

u 

to  Q. 

>-  o. 
IV  iv 

l. 

I-  >* 
IV  X 


L. 

4) 

X 


c 

o 

Q- 

3 


tO 

-O 

c 

4) 

Q. 

4> 

T3 


tO 

03 

C 


5 

C 4) 

Li 

TJ 

4- 

*4 

X L. 

c 

u 

o 

*4 

Q. 

w 

o 

Li 

o 

41 

U 

c 

O IV 

Z4 

i 

4) 

IV 

X 

&. 

4) 

O 

i/i 

f 

c 

03 

c c 

i. 

03  IV 

*4) 

"O 

4) 

X 

E 

Xi 

< 

o 

— O 

c x 

•-» 

(0 

3 

$ 

■ 

L.  — 

•— 

Li 

IV 

C 

« 

4)  4) 

4)  +* 

> 

03  C 

4) 

c 

Vj 

X U 

Li  3 

a 

m — 

3 

u 

*o 

• 

o 

c 

c 

Li  •— 

o 

E 

E 

O 

c 

*o 

Ui 

4 

TJ  4) 

iv  O 

o 

•—  4) 

c 

IV 

4) 

3 

10 

Lrf 

§ 

4)  U 
Li  10 

o 10 

40  4) 

10 

o 

(0 

M-  IV 

4- 

O 

Li 

(0 

3 

IV  4) 

fig 

w 

O 4) 

X 

10 

i 

'i 

ui 

41 

ac 

L L. 

3 O 

i X 

Li 

o 

L. 

4)  O 

4) 

10 

Li 

o 

u 

to 

(J 

Li  3 

E 03 

4) 

(0  C 

3 

2 

•fc 

UI 

H- 

> 

2 

IV  — 

(/)  u. 

2 X 

a 

«o 

3 — 

-fc 

•K 

* 

«*c 

■*? 

-38- 


J 


to 

LU 

O 

pa 

LA 

O 

CM 

UJ 

-4 

IA 

LA 

0C 

1 

i 

1 

1 

UJ 

<n 

CM 

vO 

*— 

U. 

UJ 

0£ 

CO 

-4* 

-4* 

LA 

E 

3. 

IA 

U — 

4> 

to  v 


CO  < E O 

uj  o 3 o 

CD  — CM  CM 


1- 

c 

c 

C 

o 

c 

o 

O 

4-* 

4 -» 

4J 

h- 

3 

3 

3 

CO 

CO 

-O 

XI 

4) 

jQ 

O 

»- 

•— 

•— 

N 

•• 

z 

z 

UJ 

V. 

i. 

•— 

u 

CD 

UJ 

-J 

*-» 

■M 

to 

4-» 

< 

X 

CO 

to  >* 

to  >* 

to 

UJ 

— 

— 4J 

— W 

41 

•— 

Q 

at 

CO 

T3  — 

■a  — 

ai 

“O 

CO 

O 

o 

4) 

-J 

CO 

o 

a>  o 

4)  O 

u 

4) 

< 

CJ 

2 

a. 

N — 

— 4) 

N — 

— 4) 

4) 

> 

N 

r 

CO  > 

CO  > 

< 

CO 

»- 

CL 

o 

C7> 

c 

• 

■a 

•— 

-4 

4> 

L. 

4-» 

4) 

UJ 

to 

4-> 

— J 

mm 

LJ 

CO 

4) 

IV 

< 

u. 

o 

L. 

u 

CO 

4) 

O 

>* 

— >* 

o 

Q. 

u 

>* 

CL  U 

O 

IV 

UJ 

-C 

CL  W 

c 

o 

O' 

CL 

IV 

& i 

o 

to 

o 

3 

O) 

mm 

u 

i — 

to 

u 

c 

z 

o> 

i.  a 

3 

LJ 

< 

X 

o 

4)  O 

u- 

o 

o 

to  «— 

4- 

4) 

4> 

UJ 

o 

<0  4) 

•— » 

Ql 

z 

-1  > 

O 

t/> 

X 

rr 


significant  result  of  examining  this  list  is  that  of  the  six  properties  of 
particles  listed  as  "needed"  in  2.2.,  only  two  (velocity  and  size  distribu- 
tion) are  measurable  by  the  optical  scattering  or  interference  techniques 
under  development.  The  fact  is  that  optical  scattering  techniques  can  give 
only  size  and  velocity  information  since  the  scattering  of  an  optical  beam 
is  influenced  only  by  the  scattering  cross  section  (size)  and  the  Doppler 
frequency  shift  (velocity).  Other  principles  of  physics  must  be  employed 
to  obtain  information  on  composition,  mass,  and  physical  state.  Moreover, 
the  size  range  and  the  velocities  expected  in  rocket  exhausts  are  not  fully 
addressed  by  the  optical  techniques.  Thus,  research  is  required  to  develop 
the  particle  diagnostics  needed  for  rocket  combustion  flows.  The  research 
needed  is  a re-examination  of  the  optical  scattering  techniques  for 
geometrical  size  and  velocity  and  an  investigation  of  new  approaches  for 
mass,  composition,  physical  state,  and  thermal  properties. 

3.2.2  Assessment  of  Optical  Techniques  for  Particle  Diagnostics. 

Many  investigators  are  currently  engaged  in  development  of  laser  scattering 
diagnostics  for  particle  measurements,  and  it  would  be  presumptious  of  the 
writer  to  attempt  to  survey  all  this  work.  It  will  suffice  here  to  make 
some  appropriate  comments  regarding  application  of  the  various  techniques 
to  rocket  flow  problems.  The  basic  principles  of  the  methods  are  again 
outlined  in  Appendix  B. 

3.2.2. I Holography  (Appendix  B-8) . 


In  the  strictest  sense,  holography  does  not  meet  the  criteria  laid  out  in 
Section  1 for  a diagnostic  — it  is  a direct  measurement.  But,  because  it 
is  a very  useful  in  situ  tool  with  wide  application  to  combustion  measure- 
ments, discussion  of  the  technique  is  included  here.  What  is  listed  in 

Table  4 as  limitations  on  size  and  velocity  represent  approximate  limits  on 
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optical  resolution  via  the  Rayleigh  criteria,  ' practical  limits  on 
particle  identification  for  a 1 ysec  laser  pulse  separation,  respectively. 
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The  laser  pulse  width  also  has  the  effect  of  blurring  the  particle;  for  a 

2p  particle  moving  at  200  m/sec  a pulse  width  less  than  10  nanosec  is 

required  to  "freeze"  the  particle.  These  theoretical  limits  have  been 

approached  recently  in  a system  being  constructed  for  application  to 
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combustion  bomb  type  experiments.  Extension  of  holography  to  the 
velocities  (p~3000  m/sec)  expected  in  the  near  field  of  rocket  plumes  and 
to  particles  having  dimension  less  than  2 pm  is  not  likely.  However, 
application  to  combustion  studies  and  to  the  secondary  smoke  region  (see 
Figure  1)  of  tactical  missile  plumes  can  provide  very  useful  information. 


3. 2. 2. 2 Laser  Doppler  Velocimetry  (LDV)  (Appendix  B-9) 


In  the  most  recent  approaches  to  LDV  measurements,  velocity  and  particle 

size  are  obtained  by  interpreting  the  fine  structure  of  the  scattered 

heterodyne  signal  from  the  interference  pattern  of  a probe  volume  formed 

by  two  intersecting  laser  beams  (see  Appendix  B-9).  The  velocity  is 

determined  by  the  frequency  of  the  scattering  signal  and  the  particle  size 
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by  the  "visibility"  or  shape  of  the  wave  pattern.  Two  serious  limita- 

tions arise  in  application  of  LDV  to  rocket  exhausts.  First,  electronic 
signal  processors  are  limited  to  frequencies  less  than  about  100  MHz,  which, 
in  turn,  places  an  upper  limit  of  about  1.5  km/sec  on  velocity**^  (see 
Appendix  B~9) . Nozzle  exit  and  near  field  rocket  exhaust  velocities  range 
up  to  about  3 km/sec.  (Frequency  detection  to  much  higher  frequencies  is 
possible  using  Fabry-Perot  interferometers,  but  the  practical  difficulties 
associated  with  use  of  these  devices  and  inherent  resolution  limits  negates 
their  use.  J)  Second,  the  large  particle  loading  of  solid  propellant  motor 
exhausts  makes  penetration  of  the  flow  field  difficult,  thus  seriously 
reducing  the  number  of  usable  scattering  events  for  particle  sizing  and 
velocity  measurement.  The  ultimate  lower  limit  of  measurable  particle  size 
using  this  technique  is  reached  when  the  "visibility"  (integral  of  the  pulse 
shape)  does  not  change  appreciably  with  decrease  of  particle  size.  A 

practical  lower  limit  of  about  the  wavelength  of  the  illuminating  laser 
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results  from  this  criteria.  The  upper  limit  on  size  is  governed  by  the  double 


value  of  the  Fraunhofer  diffraction  function;  for  a typical  system,  the 
upper  limit  is  about  15  pm.  Other  limitations  of  the  technique  are  due 
to  distortion  of  the  signal  pattern  by  nonhomogeneous  particle  velocity 
vectors  and  multiple  scattering  events.  These  limitations  lead  to  the 
conclusion  that  LDV  measurements  in  solid  propellant  combustion  applications 
are  restricted  to  phenomenological  laboratory  studies  and  far  field  smoke 
studies. 

3 . 2 . 2 . 3 Diffusion  Correlated  Spectroscopy  (Appendix  B-10). 

46-l«8 

A technique  that  has  been  recognized  by  only  a few  investigators  makes 

use  of  the  influence  on  the  homodyne  scattering  spectrum  from  an  LOV  system 

of  the  correlation  of  positions  of  particles  undergoing  Brownian  motion  in 

a gas.  The  technique  has  been  identified  as  "Laser  Doppler  Spectroscopy" 
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by  Hinds  and  Reist,  but  use  of  this  name  is  misleading  since  it  is  not  the 
Doppler  effect  which  determines  the  line  shape  but  the  fluctuation  in  the 
scattered  signal  received  on  a detector  caused  by  changes  in  the  phase 
factor  with  time  as  a particle  moves  via  its  Brownian  motion.  The  combina- 
tion of  all  the  scattering  signals  received  by  the  detector  from  the 
illuminated  field  results  in  a cancellation  from  noncorrelated,  independent 
particles,  but  signals  from  the  same  particle  are  correlated  and  thus  give 
a finite  correlation  function.  If  the  illuminated  field  is  uniform  and  use 
is  made  of  the  probability  of  finding  a particle  whose  motion  is  governed 
by  diffusion  at  some  position,  then  the  resulting  homodyne  power  spectrum 
will  be  given  by  a Lorentzian  line  shape.  If  the  illumination  field  is  not 
uniform,  then  the  signal  will  be  convoluted  with  the  intensity  variation 
across  the  field;  for  a Gaussian  distribution  of  intensity  (i.e.,  the  usual 
distribution  for  a single  mode  gaseous  laser  beam),  the  result  of  the 
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convolution  will  be  a Voight  profile.  The  result  of  a measurement  of  the 
line  profile  halfwidth  in  a monodjsperse  aerosol  medium  is  to  determine  the 
diffusion  coefficient,  which  is  then  uniquely  related  to  particle  size. 

For  a polydisperse  aerosol,  the  halfwidth  measurement  gives  an  "average" 
particle  size.  It  appears  that  the  theory  can  be  manipulated  to  give  a 
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distribution  of  particle  sizes  from  multiple  measurements  at  several  angles, 
but  no  serious  effort  has  yet  been  made  to  accomplish  this. 

The  potential  of  this  technique  has  only  been  scratched,  with  measurements 

i|7 

reported  only  for  laboratory  flames.  The  lack  of  application  of  the 
technique  probably  stems  from  a lack  of  understanding  of  the  principles 
involved  since,  by  the  above  discussion  and  the  literature  cited,  it  is 
clear  that  the  mathematical  physics  leading  to  the  final  result  is  quite 
complex.  Nevertheless,  the  experiment  is  very  straightforward  (see 
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Appendix  B-IO)  and  yields  much  information.  Some  work  conducted  recently 
has  contributed  to  more  confidence  in  the  method  in  that  comparison  of  the 
average  particle  size  determined  from  the  homodyne  power  spectral  halfwidth 
with  that  obtained  from  electron  microscope  photographs  of  samples  from  an 
acetylene  flame  were  in  reasonable  agreement. 

The  major  limitation  of  the  method  for  rocket  exhausts  appears  to  be 
imposed  by  the  velocity;  that  is,  a particle  must  remain  in  the  illuminated 
field  long  enough  for  the  correlation  to  be  established.  Penner,  Bernard, 
and  Chang"*®  have  estimated  this  minimum  time  to  be  about  10  ^ sec;  thus, 
for  a stream  velocity  of  1000  m/sec,  the  diameter  of  the  illumination  beam 
should  be  about  1 cm.  In  turn,  this  sort  of  requirement  necessitates  very 
high  power  lasers  in  order  to  obtain  good  signal-to-noise  ratio. 

From  the  above  discussion,  it  is  concluded  that  additional  research  on  the 
Diffusion  Correlated  Spectroscopic  Technique  is  warranted  in  order  to 
determine  its  usefulness  in  rocket  propulsion  applications. 

3.2.2.A  Mie  Angular  Scattering  (Appendix  B-ll). 

The  dependence  of  the  measured  scattering  intensity  at  a given  azimuth,  with 
respect  to  the  direction  of  an  incident  laser  beam  on  particle  size  and 
index  of  refraction  of  the  medium,  is  used  in  this  technique.  Measurements 
are  made  at  a number  of  scattering  angles  to  form  a measurement  vector. 

The  measurement  vector  can  be  equated  to  a matrix  of  angle  and  size  dependent 
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scattering  function  times  a size  vector.  The  inversion  of  the  scattering 

matrix  then  yields  a particle  size  distribution  if  the  index  of  refraction 
C|  _co 

is  known.  Presumably,  the  technique  can  also  yield  effective  indices 

of  refraction  if  some  other  measurement,  say  polarization  ratio,  is  also 
measured.  The  technique  has  not  been  thoroughly  demonstrated  with  experi- 
mental data,  but  since  it  is  based  upon  first  principle  physics,  it  should 
be  a very  powerful  technique.  Limitations  for  rocket  flows  include 

(1)  Lack  of  index  of  refraction  data  for  particles  at  the  tempera- 
tures of  interest;  bulk  material  data  is  not  adequate. 


(2)  Lack  of  sensitivity  of  the  angular  dependence  of  the  Hie 
scattering  function  on  particle  size  d in  both  the  lower 
(d<0.25A)  and  upper  (d>~10A)  size  ranges,  where  A is  the 
wavelength  of  the  incident  laser  beam. 

(3)  Only  single  scattering  data  can  be  inverted  successfully  so 
that  the  technique  is  limited  to  flows  with  particle  loadings 
sufficiently  small  that  multiple  scattering  effects  are 
negligible.  A quantitative  value  of  particle  density  which 
bounds  this  limitation  is  not  available. 


3.2.3 


Assessment  of  Gas  Phase  Diagnostics  in  Gas/Particle  Flows. 


Very  little  can  be  said  about  the  success  or  the  limitations  of  the 
spectroscopic  techniques  discussed  in  3.1  (for  clean  gas  flows)  when 
applied  to  gas/particle  flows.  Clearly,  if  the  radiation,  absorption, 
and/or  scattering  by  particles  does  not  appreciably  distort  the  emitted, 
transmitted,  and/or  scattered  radiation  from  the  gas  species,  then  the 
techniques  discussed  in  3.1  may  be  applied.  This  situation  may  exist  for 
flows  with  low  particle  loading  such  as  may  be  found  in  the  exhausts  of 
"smokeless"  type  solid  rocket  propellants,  but  experiments  to  determine 
the  extent  of  the  interference  for  such  propellants  have  not  yet  been 
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carried  out.  Some  possible  sources  of  interference  between  gas  phase 
diagnostics  and  particles  that  can  be  foreseen  are  listed  as  follows: 

(1)  Radiation  from  hot  gases  may  be  much  smaller  than  from 
particles  so  that  gas  radiation  cannot  be  detected.  (This 
situation  was  encountered  at  the  exit  of  a highly  aluminized, 
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second  stage  rocket  motor.) 

(2)  Absorption  and  scattering  by  particles  may  be  so  great  that 

no  source  radiation  is  transmitted  through  the  gas/particle 
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flow,  i.e.,  total  extinction. 

(3)  The  scattered  radiation  from  a laser  beam  may  be  multiply 
scattered  so  that  the  integrity  of  the  signal  at  the  detector 
is  lost;  i.e.,  the  detector  receives  contributions  due  to 
scattering  by  particles  all  along  the  beam  so  that  the  source 
of  the  detected  radiation  is  not  localized. 

(4)  An  incident  laser  beam  may  be  absorbed  by  the  particles, 
raising  their  temperature  and  producing  thermal  radiation 
which  contributes  to  the  background  noise  for  Raman  or 
fluorescence  detection.  (This  is  the  situation  reported  by 

i, 

Eckberth  for  Raman  measurements  in  sooty  flames.) 

The  above  list  of  possible  problems  is  not  exhaustive,  but  it  seems  clear 
that  not  enough  is  currently  known  about  how  to  avoid  or  account  for  the 
interferences  of  signal  detection  due  to  particles  to  make  gas  diagnostic 
measurements  in  highly  loaded  gas/particle  flows  with  any  degree  of  confi- 
dence. Experimental  investigation  in  simulated  gas/particle  flows  is 
warranted.  Such  investigation  should  at  least  include  determination  of 
particle  loading  and  size  bounds  under  which  the  various  gas  diagnostics 
will  not  be  affected  by  interference  caused  by  particle  effects. 
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ADVANCED  TECHNIQUES. 


r 


4. 

This  section  is  devoted  to  examination  of  techniques  having  the  inherent 
capability  for  measurement  of  those  properties  of  combustion  products  for 
which  no  measurement  techniques  currently  exist,  such  as  particle  composi- 
tion, mass,  state,  and  temperature.  Since  optical  scattering  can  at  the 
most  give  particle  size,  velocity,  and  index  of  refraction  information,  it 
is  concluded  that  other  physical  principles  must  be  exploited  to  provide 
the  other  particle  properties.  For  gas  phase  properties,  the  spectroscopic 
methods  discussed  in  3.1  contain  the  necessary  information,  but  distortions 
of  the  spectral  data  due  to  multiple  scattering  and  particle  radiation  and 
absorption  may  destroy  the  integrity  of  the  methods.  Approaches  to  both 
the  problems  of  particle  and  gas  phase  diagnostics  will  be  discussed  in  the 
following  paragraphs. 

4. 1 Particle  Composition  and  Mass,  X-Ray  Methods. 

The  specificity  of  X-ray  spectroscopy  in  identifying  and  counting  atoms  in 
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the  solid  state  is  immediately  obvious  from  the  inherent  X-ray  physics. 

An  incident  X-ray  beam  at  energy  e.  into  a field  of  particles  containing 
atoms  having  an  absorption  edge  with  energy  ej  less  than  e.  will  be 
absorbed  by  those  atoms  and  some  of  them  will  re-emit  (fluoresce)  at  an 
energy  e.  just  less  than  the  absorption  edge  energy.  The  energy  ej  or  ej 
(or  wavelengths  Aj  and  X^)  identifies  the  atom,  and  the  number  of  X-ray 
photons  emitted  (or  absorbed)  is  a measure  of  the  number  density  of  j-type 
atoms,  n^  (see  Figure  4).  Thus,  in  principle,  a method  has  been  identified 
for  particle  composition  and  mass  measurement.  However,  the  feasibility 
of  actually  utilizing  such  a method  requires  considerable  investigation. 

The  results  of  a paper  study  on  the  feasibility  of  the  use  of  X-ray  spectro- 
scopy for  gas/particle  flow  diagnostics  is  reported  elsewhere.^ 

Three  possible  X-ray  measurement  approaches  are: 


/ 

/ 
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Technique 


Possible  Measurements 


1. 

Absorption 

Average  j atom  density 

2. 

Average  fluorescence 

Average  j atom  density 

3. 

Particle  fluorescence 

Number  of  j atoms  in  a 

particle 

particle  diameter,  and 

velocity 

The  envisioned  techniques  and  the  associated  physics  are  outlined  very 
briefly  in  the  following  paragraphs. 


1* . 1 . 1 Absorption  - Average  Atomic  Density. 


When  a beam  of  X-ray  photons  of  energy  e.  is  passed  through  a particle- 
laden gas  flow,  the  attenuation  is  equal  to  the  integral  along  the  path  of 
the  sum  of  the  products  of  the  number  density  times  the  absorption  cross 
section  for  each  atomic  species,  j,  present  in  both  gas  and  particles. 
Absorption  cross  section  data  is  plentiful  for  the  atomic  species  of 
interest  in  combustion  products.  Calculations  of  the  transmittance  of 
radiation  near  the  wavelength  of  the  absorption  limit  of  aluminum  through 
the  exhaust  of  a typical  highly  aluminized  propellant  rocket  (5  cm  dia 
nozzle  exit)  are  reported  in  Reference  55.  The  following  conclusions  were 
drawn  from  the  results  of  that  study. 


(1) 


The  value  of  the  transmittance  (r  * 0.25  to  0.9  over  the 
1.4  to  10  ke  v energy  range)  is  in  a measurable  range,  and 
the  contribution  from  A i or  Cl  atoms  is  sufficient  to  measure 
the  total  number  of  A£  and  Cl  atoms  in  the  line  of  sight. 

The  contribution  from  C,  N,  and  0 elements  to  the  absorption 
must  be  known,  however,  in  order  to  uniquely  determine  the 
A£  and  Cl  content;  or  sufficient  measurements  at  different 
wavelengths  must  be  made  to  invert  the  matrix  and  thus 
determine  the  number  density  of  each  contributing  atom. 
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(2)  The  contribution  to  the  absorptance  by  0 and  N atoms  makes 
it  necessary  that  the  X-ray  path  external  to  the  rocket  flow 
be  purged  by  some  low  absorption  gas,  say  helium. 

(3)  The  Beer's  law  form  of  the  equation  for  the  X-ray  absorption 

is  amenable  to  Abe)  inversion  so  that  local,  average  density 

can  be  determined  for  axisymmetric  flows  from  the  measured 

transmittance  profiles.  The  three-dimensional  inversion 

techniques  for  X-ray  absorption  as  used  by  the  medical 
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profession  also  appear  to  be  adaptable  for  nonaxisymmetric 
cases. 

The  above  conclusions  drawn  from  Reference  51  provide  sufficient  justification 
to  conduct  an  experimental  demonstration  study  of  the  X-ray  absorption 
techniques  with  consideration  of  available  equipment  and  definition  of 
rocket  test  scenarios.  It  should  be  emphasized  that  this  technique  can 
only  yield  an  average  atomic  density,  not  direct  information  about  particles. 
However,  if  the  A l is  contained  in  the  A^O^  particles  as  believed  by  most 
combustion  experts,  then  a measure  of  the  local  average  density  distribution 
can  yield  much  information  about  the  validity  of  the  gas-dynamic  treatment 
of  two-phase  flows  and  the  assumed  particle  size  distributions  which  are 
input  to  the  models. 

A . 1 . 2 Time  Averaged  Fluorescence  - Space  Averaged  Atomic  Density. 

Many  of  the  problems  associated  with  the  X-ray  absorption  measurement  (i.e., 
interferences  from  other  elements,  inversion  to  obtain  local  values,  etc.) 
are  avoided  if  the  characteristic  fluorescent  radiation  at  the  Ka  energy  of 
the  specie  of  interest  can  be  measured.  An  array  of  detectors  tuned  to  the 
energy  e^  can  then  be  used  to  obtain  local  values  of  the  atomic  number 
density.  The  main  problems  with  the  fluorescence  measurements  are  the 
determination  of  the  local  beam  intensity  and  the  absorption  in  the  path 
between  the  X-ray  beam  and  the  detectors.  The  rates  of  Ail  and  Ct  Ka  photon 


emission  at  the  axis  of  a typical  rocket  and  the  absorption  along  the  path 

to  the  detector  were  calculated  and  are  given  in  Reference  55.  It  was 

shown  that  the  fluorescence  measurement  can  be  accomplished  provided  that 

1 2 2 

X-ray  sources  having  a flux  on  the  order  of  10  photons/sec  cm  are 
provided.  Again,  the  theoretical  calculations  show  sufficient  merit  that 
laboratory  demonstration  experiments  are  justified. 

i* . 1 . 3 Particle  Fluorescence  - Number  of  Atoms  in  Particle,  Velocity,  Size. 

The  maximum  yield  of  information  about  particles  using  X-ray  fluorescence 
would  come  from  the  pulse  of  photons  produced  as  a particle  moves  through 
the  X-ray  beam.  Conceptually,  the  fluorescent  event  could  yield  information 
on  the  kind  and  number  of  atoms  in  the  particle,  the  mass,  the  velocity,  and 
the  geometrical  dimension  (see  Figure  5).  However,  the  same  sort  of  problems 
that  are  encountered  in  optical  scattering  (multiple  scattering,  finite 
illuminated  volume,  etc.)  are  encountered  here,  as  well  as  the  same  absorption 
problem  that  accompanies  the  interpretation  of  average  X-ray  fluorescence 
measurements  (4.1.2).  The  total  possible  number  of  A£.  Ka  photons  emitted  by 
a Ip  AJZ^O^  particle  traveling  at  1000  m/sec  during  the  crossing  of  a 1 mm 
wide  X-ray  beam  (~1  m/sec)  of  X-ray  photons  at  an  energy  just  larger  than 
the  absorption  edge  is  estimated  to  be  10  times  the  incident  X-ray  photon 
flux  rate.^  Some  30  percent  of  these  will  be  absorbed  in  the  path  to  the 
detector  in  the  typical  tactical  rocket  exhaust  flow.  The  questions  to  be 
answered  by  a feasibility  study  are  then,  "Does  there  exist  a source  of 
X-rays  of  sufficient  intensity  to  produce  about  100  photons  from  the  particle 
in  the  one  more  time  interval?"  "What  are  the  detector  requirements  to  record 
a pulse  shape  (or,  less  desirably,  the  integral  over  the  pulse)?";  "Does  such 
a detector  exist,  and,  if  not,  is  such  a detector  possible?"  Again,  the 
potential  of  this  method  is  sufficient  great  that  laboratory  investigation 
is  justified. 


-50- 


At  = d/v 

SELECTED  SINGLE  EVENTS 


Figure  5.  Scattering  Event  Chronology  for  a Particle 
Passing  Through  an  X-Ray  Beam. 


It. 2 Particle  Temperature  - Radiation/Absorption  Properties. 

The  usual  approach  to  the  measurement  of  the  temperature  of  a particle 

cloud  that  is  not  optically  thick  is  to  determine  the  spectrum,  look  up 

an  emissivity  value  for  the  particle  material,  and  using  the  measured 

radiation  and  the  Planck  function  at  two  or  more  wavelengths,  calculate  a 

temperature.  This  is  two  color  pyrometry;  or,  for  optically  dense  clouds, 

an  optical  pyrometer  will  do  the  Job  at  one  wavelength.  Still  another 

approach  which  can  be  used  if  the  emissivity  is  well  behaved  is  the  use  of 

the  Wein  displacement  law,  AT*  2898pm  K.  Two  problems  arise  with  these 

m 

approaches: 

(1)  Emissivity  data  for  bulk  material  is  not  the  same  as  for 
particles,^  and 

(2)  the  Planck  function  (and  Wein  displacement  law)  must  be 
modified  when  the  particle  dimension  is  close  to  the  value 

rO 

of  the  wavelength  of  the  radiation . 

The  first  problem  is  unsolved  due  to  lack  of  experimental  data.  The  second 
problem  has  recently  been  recognized,  and  a revised  formulation  for  equili- 
brium radiation  from  black  particles  has  been  derived  but  not  experimentally 
verified.^ 

Thus,  two  very  important  areas  of  research  are  identified,  and  measurement 
of  particle  temperature  must  await  solution  of  these  problems. 

Stated  very  briefly,  the  revision  of  the  Planck  Black  Body  Law  for  particles 

of  dimension  close  to  the  wavelength  of  interest  can  be  expressed  (to  first 

order)  in  the  form  of  a correction  term.  For  spherical  particles,  the 

revised  function  Nbb  (A,  T,  d)  is  dependent  upon  wavelength  A,  temperature 

T,  and  the  particle  diameter  d,  and  is  related  to  the  Planck  function  N_ 
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(A,  T)  through  the  expression  ' 
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where  A cannot  be  greater  thanV5~7rd.  Figure  6 gives  a typical  result  for 
1 000° K particles  over  a range  of  particle  sizes.  The  result  is  quite 
striking  and  obviously  affects  the  determination  of  temperature  from  the 
spectral  shape  of  the  radiation  from  particles  or  the  Wein  displacement 
law.  The  portion  of  the  spectrum  at  wavelengths  less  than  about  the  diameter 
of  the  particle  are  not  appreciably  altered,  however,  and  the  temperature  of 
the  particles  may  then  still  be  obtained  by  the  Planck  function.  Experimental 
verification  of  the  revised  radiation  law  is  sorely  needed. 

The  work  by  Dowling  and  Randall"^  addresses  the  measurement  of  small  particle 
radiation,  but  results  have  not  yet  been  forthcoming.  In  their  experiment, 
the  radiation  from  a single  particle  trapped  in  the  focal  region  of  two 
laser  beams  is  sought.  The  validity  of  the  revised  radiation  expressions 
can  be  tested  by  the  dependence  of  the  radiation  spectrum  on  different  size 
particles  using  their  experiment.  A prerequisite  to  measurement  of  temper- 
ature of  particles  is  resolution  of  the  particle  radiation  law  and  the 
emissivity  problems.  If  the  experiment  by  Randall,  et.  al.,  does  not  yield 
information,  then  experiments  using  particle  clouds  with  all  the  incumbent 
problems  of  particle  coagulation  are  necessary  in  order  to  develop  an  in 
situ  measurement  technique  for  particle  temperature. 


**. 3 Particle  State  - Liquid  or  Solid. 

Evidence  exists  to  show  that  AJt-0-  particles  are  in  the  liquid  state  over 
at  least  a part  of  their  trajectory  in  a rocket  motor.  This  conclusion 
can  be  drawn  in  two  ways.  First,  the  temperature  within  the  motor  and 
through  the  initial  expansion  in  the  nozzle  is  much  greater  than  the  melting 
temperature  of  (i.e.,  about  2200°K).  Second,  experimental  evidence 

in  the  residue  on  nozzle  walls  and  in  the  construction  of  the  larger 
particles  that  are  collected  after  a firing  points  to  a prior  liquid  state. 


i 


1 


\ 
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But  the  information  sought  for  use  in  performance  models  is  the  location 
in  the  nozzle  where  the  state  change  takes  place.  If  a diagnostic  were 
available,  the  information  could  presumably  be  obtained  by  a series  of 
tests  using  different  nozzle  expansion  ratios. 

One  principle  that  can  be  cited  as  a possible  diagnostic  is  the  expected 
change  of  refractive  index  as  the  particle  changes  state,  which  should  be 
evidenced  in  the  spectral  radiation  characteristics  through  the  apparent 
emissivity  or  in  the  Mie  scattering  characteristics.  In  either  case, 
insufficient  information  exists  to  propose  a diagnostic,  and,  again, 
research  is  needed  on  the  radiation  characteristics  of  particulate  matter. 

A.  4 Gas  Properties  - Spectroscopic  Methods. 

All  the  techniques  discussed  in  3*1  have  the  potential  of  being  successful 
in  gas/particle  flows  under  particle  loading  conditions  in  which  the  effect 
of  particles  on  the  gas  induced  radiation  are  minimal.  It  remains  to 
determine  the  conditions  under  which  a given  technique  can  be  used,  by  both 
modeling  the  process  and  conducting  laboratory  experiments.  A suggested 
research  plan  is  then  very  straightforward: 

(1)  Model  the  radiation-particle  interaction  to  predict  the 
effects  of  scattering  and  particle  radiation  and  absorption 
on  the  emitted,  transmitted,  or  scattered  radiation  of 
interest. 

(2)  Under  controlled  conditions  of  both  gas  and  particle  pro- 
perties, conduct  experiments  to  provide  confidence  in  the 
model  predictions  and  to  obtain  experience  with  the  inter- 
ferences on  signals  resulting  from  particles  in  the  flow. 

(3)  From  the  results  of  steps  1 and  2,  determine  optimum 
diagnostic  techniques  to  use  for  different  types  of  partic- 
ulates and  different  size  and  density  ranges. 
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It  is  recognized  that  the  recommended  approach  presented  above  is  not  a 
very  glamorous  one,  but  no  devious  schemes  to  aviod  particle  interferences 
have  so  far  come  to  mind.  Perhaps  the  straightforward  attack  suggested 
will  serve  to  inspire  ingenuity. 

Considering  the  spectroscopic  methods  available,  the  following  specific 
comments  can  be  made. 

(1)  The  low  signal  levels  encountered  in  Raman  scattering  make 
this  technique  very  susceptible  to  background  emission. 

Mie  scattering  from  particles  is  several  orders  of  magnitude 

greater  than  Raman  scattering  from  molecules  and  presents  a 

serious  background  problem.  Also,  the  experience  of  Eckberth 
L 

et.  at.,  with  laser  heating  of  the  particles  tends  to  make 
the  technique  even  more  questionable  for  hot  gas/particle 
flows.  Eckberth  concluded  that  resonant  enhanced  methods 
such  as  Coherent  Antistokes  Raman  Scattering  (CARS)  in  which 
much  larger  signal  levels  can  be  expected,  should  be  used. 
However,  the  nonlinear  dependence  on  both  beam  intensity  and 
molecular  density  have  discouraged  this  writer  from  serious 
consideration  of  CARS  for  solid  rocket  propulsion  applications 
(see  also  Reference  61). 

(2)  Methods  which  require  absorption  measurement  are  susceptible 
to  error  in  that  the  particle  contribution  to  absorption  and 
scattering  must  be  taken  into  account.  Moreover,  beam  defini- 
tion becomes  problematical  for  high  particle  laden  flows. 

(3)  Saturated  fluorescence  methods  appear  promising,  with  the 
main  problem  being  that  of  localizing  the  observation  point 
to  reduce  effects  of  multiple  scattering. 
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(1»)  Most  solid  propellants  contain  traces  of  sodium  and  potassium 
impurities.  Because  of  the  low  ionization  potential  of  these 
alkali  atoms,  ionization  and  emission  of  radiation  from  upper 
excited  electronic  status  occur.  If  the  spectra  from  the 
upper  excited  states  of  alkali  metals  is  observable,  the 
temperature,  the  electron  concentration,  and  the  alkali  atom 
density  can  be  determined  from  emission  measurements. 


5. 


CONCLUSIONS. 


The  study  reported  herein  has  considered  the  needs  for  measurement  of 
properties  of  combustion  products  In  rocket  motors  and  exhaust  plumes. 
The  conclusions  reached  during  the  study  are  listed  as  follows: 


(1)  The  needed  measurements  for  rocket  combustion  serve  mostly 
as  inputs  and  validation  data  for  propulsion  system  perfor- 
mance evaluation,  internal  models  of  combustion  and  nozzle 
expansion  processes,  and  prediction  models  for  plume  observ- 
ables (e.g.,  IR  radiation,  plume  visibility,  radar  cross 
sections).  The  desired  measurements  include  gas  temperature 
and  major  species  concentration,  and  the  temperature,  size, 
mass,  composition,  state  and  velocity  of  particles  in  flows 
ranging  from  3500°K  - 1*0  atm  rocket  chamber  conditions  to 
the  low  density  high  altitude  plumes  at  ambient  temperature 
and  densities. 


(2)  For  the  case  of  purely  gaseous  flows  from  modern  liquid 
propellant  rockets,  the  currently  existing  spectroscopic 
techniques  of  measurement  (including  high  and  low  resolution 
spectroscopy,  resonance  absorption,  fluorescence,  electron 
beam  fluorescence,  and  laser  Raman  scattering)  were  found  to 
be  adequate,  in  principle,  to  meet  rocket  propulsion  system 
needs.  Refinements  in  the  form  of  instrumentation  improve- 
ment, data  treatment,  and  sensitivity  and  uncertainty  analysis 
will  be  required  on  a continuous  basis  to  meet  increased 
accuracies  required  in  the  validation  of  models,  and  the  need 
for  more  accurate  fundamental  data  on  molecular  parameter, 
cross  section,  etc.,  is  ever  present. 

(3)  For  the  case  of  the  combined  gas/particle  flows  from  solid 
propellant  combustion,  current  techniques  for  measurement 
of  particle  properties  were  found  to  be  inadequate,  and 
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techniques  for  measurement  of  gas  properties  in  the  presence 
of  particles  were  found  to  be  largely  untried.  Optical 
techniques  for  particle  diagnostics  address  only  the  particle 
size  and  velocity  problems,  and  other  principles  must  be 
employed  to  measure  composition,  mass,  state  and  temperature. 
Gas  phase  methods  suffer  from  interferences  caused  by  back- 
ground emission  from  particles,  multiple  scattering  which 
destroys  the  geometrical  integrity  of  the  data,  and  severe 
attenuation  and  dispersion  of  source  beams. 

(4)  X-ray  spectroscopy  appears  promising  for  measurement  of 
particle  composition  and  mass  properties.  Absorption  and/or 
average  fluorescence  are  feasible  for  measurement  of  A£  or 
C£  local  atomic  densities  and  possibly  can  provide  measure- 
ment of  0,  N,  and  C densities.  From  these  measurements  and 
a suitable  model,  particle  mass  distributions  can  be  inferred. 
Individual  particle  fluorescent  scattering  events  are  a 
possibility  for  direct  measurement  of  particle  atom  content 
and  velocity.  Further  feasibility  analysis  on  all  X-ray 
methods  to  define  instrumentation  requirements  and  data 
analysis  is  required,  to  be  followed  by  proof  of  principle 
experiments. 

(5)  Only  emitted  or  absorbed  spectrally-resolved  radiation  can  be 
used  to  measure  the  temperature  of  particles.  The  emissivity 
and  radiation  law  for  particles  is  dependent  upon  the  particle 
size  when  the  size  is  near  that  of  the  wavelength  of  the 
radiation.  Expressions  derived  for  "black  particle"  radiation 
based  upon  calculation  of  the  eigenstate  density  for  sma’l 
particles  and  employing  Bose-Einstein  quantum  statistics  must 
be  validated  by  experiment,  and  particle  emissivities  for  the 
materials  of  interest  must  be  measured  before  further  progress 
can  be  made  toward  particle  temperature  measurement. 
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(6) 


A systematic  modeling  and  experimental  program  is  required 
to  investigate  spectroscopic  techniques  for  gas  property 
measurement  in  the  presence  of  particles.  The  methods  of 
fluorescent  scattering  and  alkali  metal  emission  appear  to 
offer  the  best  possibilities  and  should  be  investigated 
first. 
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APPENDIX  A.  EQUILIBRIUM  PROPERTIES  OF  COMBUSTION  PRODUCTS  FOR 
SEVERAL  TYPICAL  ROCKET  PROPELLANT  COMBINATIONS. 


This  appendix  contains  results  of  one-dimensional,  equilibrium  calculations 
of  rocket  motor  performance  and  species  concentrations  for  propellants  that 
are  typical  for  propulsion  systems  for  which  measurements  may  be  required. 
The  computer  code  used  was  an  unpublished  code  used  by  the  AFRPL  and  is 
similar  to  that  presented  in  CPIA  publication  #246,  "JANNAF  Rocket  Engine 
Performance  Prediction  and  Evaluation,"  April  1975. 


TABLE  A-1.  02/H2. 


Expansion  Ratio 

Chamber 

6: 1 

25:1 

40:1 

1 (Sec) 

sp 

0.0 

407.9 

446.3 

455.0 

Temperature  (°K) 

3107 

1678 

1109 

959 

Species 

Hole  Fractions 

H 

0.024 

- 

- 

- 

HO 

0.011 

- 

- 

- 

H2 

0.419 

0.433 

0.433 

0.433 

H2° 

0.545 

0.567 

0.567 

0.567 

TABLE  A-2. 

. lo2/rp-i. 

Expansion  Ratio 

Chamber 

6:1 

25=1 

40:1 

1 (Sec) 

310.2 

342.7 

350.5 

sp 

Temperature  (°K) 

3484 

2170 

1541 

1378 

Species 

Mole  Fractions 

CO 

0.374 

0.355 

0.321 

0.304 

CO 

0. 109 

0.151 

0.187 

0.204 

H 

0.035 

- 

- 

- 

HO 

0.038 

0.002 

- 

- 

H2 

0.130 

0.154 

0.190 

0.207 

h2o 

0.302 

0.337 

0.302 

0.285 

°2 

0.006 

- 

- 

- 

0 

0.006 

- 

- 

- 
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IABLE a- 3.  n2o1|/mmh. 


Expansion  Ratio 

Chamber 

6: 1 

25:1 

40:1 

i (Sec) 

sp 

- 

273.3 

328.5 

Temperature  (°K) 

3241 

1816 

1241 

Species 

Mole  Fractions 

CO 

0.110 

0.098 

0.071 

0.067 

C02 

0.646 

0.071 

0.099 

0.107 

H 

0.019 

- 

- 

- 

OH 

0.022 

- 

- 

- 

H2 

0.120 

0.134 

9.162 

0.177 

h2o 

0.355 

0.374 

0.345 

0.326 

NO 

0.013 

- 

- 

- 

N2 

0.313 

0.323 

0.323 

0.323 

°2 

0.002 

- 

- 

- 

TABLE  A-4.  ClFj/BjHg. 


Expansion  Ratio 

Chamber 

6;  1 

25:1 

40:1 

'sp  (Sec) 

295.  A 

331.1 

340.4 

Temperature  (°K) 

Species 

4320 

3028 

Mole  Fractions 

2618 

2509 

BCL 

.005 

.003 

.003 

.002 

BCLF 

.008 

.006 

.004 

.004 

bclf2 

.005 

.012 

.017 

.018 

bcl2 

.001 

- 

- 

- 

bcl2f 

- 

.001 

.002 

.002 

BF 

.092 

.076 

.070 

.068 

BF2 

.180 

.190 

.188 

.173 

bf2h 

.002 

.002 

.002 

.002 

BF3 

.031 

.080 

.122 

.137 

CL 

.160 

.100 

.066 

.057 

CLH 

.111 

.212 

.255 

.281 

CL2 

.001 

- 

- 

- 

F 

.021 

.001 

- 

- 

FH 

.287 

.271 

.239 

.228 

H 

O 

VO 

O' 

.01*3 

.032 

.028 
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TABLE  A-5.  ClFj/AZ-50. 


Expansion  Ratio 

Chamber 

6: 1 

25:1 

40:1 

‘sp  (S‘C» 

310.72 

326.96 

333.41 

Temperature  (°K) 

Species 

3527 

1815 

Mole  Fractions 

1359 

1159 

CLN 

.001 

- 

- 

- 

CF 

.001 

- 

- 

- 

CHN 

.029 

- 

- 

- 

CN 

.003 

- 

- 

- 

c2h 

.009 

- 

- 

- 

c2hcl 

.011 

- 

- 

- 

C2H2 

.008 

- 

- 

- 

c2n 

.002 

- 

- 

- 

LL 

.051 

- 

- 

- 

CLH 

.131 

• 179 

.179 

.179 

F 

.004 

- 

- 

- 

FH 

.542 

.536 

.536 

.536 

H 

.019 

- 

- 

- 

H2 

.032 

.037 

.037 

.037 

N2 

.133 

.146 

.146 

.146 

C 

.024 

. 102 

.102 

.102 

(Graphite) 


TABLE  A-6.  MINIMUM  SMOKE  SOLID 


• 

Expansion  Ratio 

Chamber 

111 

3.7:1 

1 (Sec) 

sp 

98.5 

201.3 

Temperature  (°K) 

2323 

2068 

1240 

Pressure  (Psia) 

1000 

555 

44.2 

Species 

Mole  Fractions 

CO 

0.458 

0.452 

0.402 

C02 

0.083 

0.089 

0. 138 

H2 

0.175 

0.181 

0.231 

H2° 

0.176 

0. 170 

0. 120 

N2 

0.108 

0.108 

0. 108 

TABLE  A-7.  ALUMINIZED  (1 8 PERCENT)  SOLID 


Expansion  Ratio 

Chamber 

',P  ««> 

- 

Temperature  (°K) 

3394 

Pressure  (Psla) 

1750 

Species 

CO 

0.240 

co2 

0.013 

H 

0.022 

H2 

0.294 

H2° 

0.125 

OH 

0.004 

M2 

0.077 

Cl 

0.007 

HC1 

0.135 

Alcl2 

0.002 

A.2°3 

0.073 

1:1 

7:1 

15: 1 

23:1 

106.4 

250.8 

274.4 

284.8 

3175 

2178 

1849 

1689 

1004 

40.2 

14.7 

8.6 

Mole  Fractions 

0.241  0.239 

0.236 

0.233 

0.013 

0.0176 

0.022 

0.024 

0.017 

- 

- 

- 

0.299 

0.317 

0.317 

0.320 

0.124 

0.120 

0.115 

0.112 

0.003 

- 

- 

- 

0.076 

0.079 

0.079 

0.079 

0.005 

0.001 

- 

- 

0.140 

0.151 

0.152 

0.152 

0.001 

- 

- 

- 

0.075 

0.076 

0.076 

0.076 

-72- 


APPENDIX  B.  BRIEF  DESCRIPTION  OF  EXISTING  DIAGNOSTIC  TECHNIQUES 
APPLICABLE  TO  ROCKET  COMBUSTION  AND  PLUME  MEASUREMENTS. 


Introduction 

The  purpose  of  this  appendix  is  to  present  in  the  briefest  way  the  essential 
principles,  assumptions,  hardware,  restrictions  and  limitations,  and  uncer- 
tainty estimates  for  the  techniques  that  have  been  employed  either  for 
combustion  related  diagnostics  or  in  similar  fields.  This  material  is 
included  in  the  report  for  completeness  and  to  support  statements  made  in 
the  main  text.  An  example  result  of  an  application  is  given  where  possible. 
The  examples  given  are  either  of  test-oriented  rocket  applications  or  as 
closely  related  to  the  rocket  measurement  as  possible.  The  choice  of  examples 
was  somewhat  arbitrarily  made  to  illustrate  the  status  of  equipment  develop- 
ment and  the  sophistication  of  the  techniques  according  to  published  material, 
and  the  reader  may  have  knowledge  of  better  examples.  References  were  chosen 
which,  in  the  author's  judgement,  best  describe  the  technique. 

The  principles  of  each  technique  are  expressed  mathematically  wherever 
possible,  calling  attention  to  "first  principle"  determination  of  the 
desired  properties  from  the  measured  variables.  The  equations  have  been 
condensed  as  much  as  deemed  possible,  and  the  reader  is  referred  to  the 
references  for  extensive  detial.  The  restrictions  and  weaknesses  listed 
stress  the  fundamental  physics  in  most  cases,  believing  that  hardware 
development  is  mostly  a matter  of  time  and  money. 

The  principal,  descriptive  reference  for  each  technique  is  called  out  in 
the  title  of  each  section,  with  other  supporting  references  introduced  as 
needed.  Some  references  are  repeated  here  from  the  main  text  in  order  that 
each  description  might  stand  alone.  Again,  only  the  references  deemed 
sufficient  to  describe  the  technique  are  listed,  and  no  attempt  was  made  to 
be  exhaustive. 
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APPENDIX  B-l.  IR  BAND  MODEL  EMISSION/ABSORPTION 
(REFERENCES  B-5  AND  B-6). 


Principle 


Emission  and  absorption  spectroscopy  can  be  applied  to  gas  flow  diagnostics 
most  straightforwardly  by  dealing  with  resolved  spectral  lines  (see  Appendix 
B-5).  If  the  liner  of  a band  system  cannot  be  resolved,  but  the  molecular 
species  has  a well  defined  structure,  then  line  by  line  calculations  can  be 
performed  (see  Appendix  B-2).  However,  for  many  species  (H^O  and  CO^,  for 
example)  the  structure  of  the  band  systems  have  not  been  defined,  and  the 
line  spacing  is  such  that  many  lines  may  fall  within  a practical  instrument 
bandpass.  For  such  cases  an  artifice  called  a "band  model"  is  employed 
(References  B-l  and  B-2),  in  which  the  ordinary  laws  of  radiative  transfer 
are  applied  to  a fictitious  line  by  line  model  which  is  designed  to  approxi- 
mate the  radiance  and  transmittance  of  the  actual  molecular  radiance  and 
transmittance  measured  at  a specified  temperature,  partial  pressure,  and 
gas  density.  Such  models  have  been  shown  to  simulate  the  actual  radiative 
transfer  quite  well  within  the  accuracy  of  the  empirical  data  available 
(Reference  B-2).  Since  H^O  and  CO2  are  the  principle  products  of  combustion 
in  most  practical  rockets  and  both  have  accessible  infrared  (IR)  emission/ 
absorption  bands  with  extensive  band  model  data  (References  B-3  and  B-4) , a 
diagnostic  technique  has  been  developed  around  the  radiative  transfer 
equations  using  the  band  model  theory  for  these  molecules  (References  B-5 
and  B-6) . 


The  first  element  in  the  diagnostic  technique  is  the  correct  expression  for 
the  radiative  transfer.  Considering  a path  through  a non homogeneous  gas 
stream  made  up  of  many  small  isothermal,  isoboric  zones  of  constant  composi- 
tion, the  radiance  NAX  emitted  at  the  edge  of  the  path  in  a finite  band 
increment  AX  is  given  by  (see  Figure  B-l  1) 


AX 


m 

£.  £ 

J 


i-1  NB  <TiAX>  *Ti-l  i ’ Pi  j ' Pi'  ‘ Ti  *Ti’  Pij’  Pi^  (B"’ 
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Figure  B-l  ].  Illustration  of  Radiative  Transfer  Through  Nonhomogeneous 

Radiating  and  Absorbing  Path  Made  up  of  Isothermal,  Isobaric 
Zones. 


IT 


where  NQ  (t.)  is  the  radiance  of  a black  body  at  temperature  T{  within 
the  wavelength  increment  AX; 


t . _ j j are  the  transmi ttances  from  the  Initial  edge  (i  ■ 1) 
to  the  beginning  of  the  I th  zone  (i  - 1)  and  to  the  end  of  the 
ith  zone  (i),  respectively; 


T.  is  the  temperature  in  the  ith  zone; 


P.  is  the  static  pressure  in  the  ith  zone;  and 


P.j  is  the  partial  pressure  of  the  ith  radiating  species  in  the 
ith  zone. 


The  transmittance  from  the  edge  i ■ 1 to  the  end  of  the  ith  zone  is  given  by 


In  (t.) 


h-1  (klAxM 


(B- 1 2) 


where  (k.  )..  represents  the  absorption  coefficient  of  the  ith  zone  made  up 

h A A 

of  contributions  from  all  the  spectral  lines  within  the  wavelength  internal 
AX,  and  depends  upon  the  temperature  T^,  the  static  pressure  Ph>  the  partial 
pressure  of  the  absorbing  species  Phj  , and  the  partial  pressure  of  all  other 
species  that  goes  to  make  up  Ph-  The  dependence  upon  other  species  comes 
about  through  foreign  gas  contributions  to  line  broadening.  The  complicated, 
unknown  character  of  the  absorption  coefficient  is  circumscribed  through  the 
artifice  of  the  band  model  and  associated  experimentally-determined  band 
model  parameters. 

In  the  nomenclature  of  Reference  B-3,  the  band  model  expression  for  t for  a 
homogeneous  path  is 


i In  (t)ax  - -2ir  (^)  f (x) 


(B-l  3) 
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where  y is  the  half  width  (width  at  half  maximum  intensity)  of  the  lines 
within  the  wavelength  increment  and  depends  upon  the  static 
pressure,  the  composition,  and  the  temperature  of  the  gas  within 
the  homogeneous  increment  l\ 

d is  the  spacing  between  the  lines; 


f (x)  is  a function  of  the  parameter 


. I /S/d, 

- 27  W 


( B— 1 A) 


which  depends  upon  the  form  of  the  band  model  being  used, 
including  the  line  shape  function; 

Z is  the  length  of  path  over  which  the  temperature,  composition, 
pressure,  and  partial  pressure  of  the  absorbing  species  are 
constant;  and 

S is  the  line  strength,  defined  for  single  lines  as 


,°°  . 

o KV 


(B-l  5) 


Now,  the  x.  and  T._ j of  Equation  (B-l  1)  are  made  up  of  contributions  from 
each  zone  along  some  nonhomogeneous  path  between  i = 1 and  i (or  i-l), 
while  the  transmittance  expressed  by  Equation  (B-l  3)  applies  only  to  a 
homogeneous  path.  A second  artifice  must  be  used  to  relate  the  defined, 
homogeneous  band  model  parameters  to  the  integrated  path  transmittance. 

One  approach  is  referred  to  as  the  Curt is-Godson  approximation  (References 
B-3  and  B-7),  which  states  that  "there  is  some  homogeneous  gas  with  band 
model  parameters  (S/d)  and  (y/d)  which  has  the  same  transmittance  as  the 
inhomogeneous  gas  having  individual  zone  parameters  (S/d)^  and  (y/d)^" 

The  result  of  the  approximation  is 


( B—  I 6) 
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where  L is  the  total  path  length  composed  of  the  sum  of  the  individual  zonal 
path  lengths,  Young  (Reference  B-4)  discusses  the  limitations  of  the 

Curt is-Godson  approximation  and  introduces  alternatives,  which  have  the 
effect  of  changing  the  form  of  f (x) . 

The  key  to  the  use  of  the  band  model  approach  to  radiative  transfer  is  an 
adequate  expression  f (x) , which  will  satisfy  experimental  data  over  a 
sufficiently  large  range  of  temperature  and  pressure.  Full  discussion  of 
the  current  forms  for  f (x)  under  all  conditions  of  temperature  and  pressure 
and  for  different  species  is  beyond  the  scope  of  this  brief  discussion,  und 
the  references  must  be  consulted  (see  References  B-3,  B-4,  B-7,  B-8,  and  B-9)  • 
The  recommended  function  for  use  with  H20  (2J-26  pm)  and  C02  (4. 35-4*- pm) 
is  that  reported  by  Young  (Reference  B-10)  which  is  based  upon  an  exponential- 
tailed line  strength  distribution  with  combined  collisional  and  pressure 
broadening. 

The  second  part  of  the  diagnostic  approach  using  the  band  model  radiative 
transfer  equations  is  the  necessary  inversion  procedure  to  determine  local 
values  of  temperature  and  species  partial  pressure  from  the  integrated 
radiance  and  transmittance  measurements.  The  procedure  developed  by  Brewer 
and  Limbaugh  (Reference  B-5)  for  cyl indrlcal ly  symmetric  flows  is  illustrated 
in  Figure  B-l  2.  The  symmetric  flow  is  divided  into  concentric  zones  in  which 
the  temperature,  static  pressure,  and  species  partial  pressure  are  assumed 
constant.  If  there  are  M such  zones  and  M measurements  of  transmittance  and 
radiance  are  made  as  shown  in  Figure  B-l  2,  and  the  local  composition  and 
static  pressure  are  known  a priori,  then  sufficient  information  is  available 
to  determine  zonal  temperatures  and  partial  pressures.  The  procedure  is  to 
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begin  at  the  outer  zone  where  there  are  two  unknowns  and  two  measurements 
and  the  T,  and  Pjj  can  be  determined  immediately  from  Equations  (B- 1 1)  and 
(B-1  2)  using  a weighted  path  length  which  depends  upon  the  width  of  the 
concentric  zone  spacing.  Then,  the  calculation  procedes  to  the  second  path 
in  which  now  the  zone  2 temperature  and  partial  pressure  are  the  unknowns 
and  an  iteration  is  required  to  converge  to  a solution  which  will  satisfy 
the  measured  radiance  and  transmittance.  The  procedure  is  repeated  for 
path  3,  path  I*,  etc.,  until  values  of  temperature  and  partial  pressure  are 
determined  for  all  the  concentric  zones  M,  as  shown  in  Figure  B-1  3.  This 
particular  procedure  is  sometimes  referred  to  as  the  "onion  peel"  process. 
Other  procedures  for  the  axisymmetrlc  inversion,  such  as  an  iterative  solu- 
tion of  the  Abel  integral  equation  (Reference  B-10)  may  be  equally  adequate. 
However,  the  method  of  solution  employed  by  Limbaugh  (Reference  B-U)  makes 
it  possible  to  conventionally  carry  propagation  of  measurement  uncertainty 
in  radiance  and  transmittance  through  the  mathematical  procedure  and  arrive 
at  a consistent  uncertainty  estimate  in  temperature  and  partial  pressure. 

Restrictions  and  Weaknesses 

1.  The  technique  is  limited  to  molecules  having  isolated  bands  with 
closely  spaced  lines  for  which  a band  model  has  been  adapted  and 
parameters  obtained.  Currently,  application  is  limited  to  H^O 
(2.5  ym)  and  COj  (**.%j). 

2.  The  mathematical  inversion  of  the  experimental  measurements  of  radiance 
and  transmittance  to  temperature  and  partial  pressure  involves  a number 
of  sensitivity  and  uncertainty  problems.  In  particular,  the  sensitivity 
of  inner  zone  properties  determined  from  inversion  of  the  experimental 
data  to  radiance  and  transmittance  is  small  compared  to  outer  zone 
properties. 

3.  The  validity  of  the  inversion  depends  upon  symmetry  of  the  source,  so 


that  applications  are  limited. 


1 
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N - RADIANCE 
t - TRANSMITTANCE 


T - TEMPERATURE 
P - SPECIES  PARTIAL  PRESSURE 


ITERATION 


Figure  B-l  3-  Computational  Flow  Chart  for  AEDC  IR 
Emission/Absorption  Inversion  Method. 
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4.  The  accuracy  of  the  band  model  parameters  is  not  well  established  so 
that  an  absolute  uncertainty  estimate  is  not  possible. 

Implementation 

A successful  measurement  scheme  Is  reported  In  Reference  B-12,  consisting  of 

Source:  calibrated  1.27  cm  aperture  black  body  chopped  at 

400  Hz  focused  by  Ca  lens  at  center  of  flow  field. 

Receiver:  matching  Ca  F2  lens 

chopper  at  333  Hz 

interference  filters 

-H20:  54  cm  ' half  width  centered  at  4005  cm  * 

C02:  36  cm  1 half  width  centered  at  2278  cm  ' 

Detector:  In  Sb 

Electronics:  two  lock-in  amplifiers,  one  synchronized  with 

source  chopper  and  one  with  the  receiver  chopper 


Data  recorded  on  digital  tape 

Examp  1 e 

T and  [CO^]  profiles  were  obtained  from  the  exit  plane  of  a Transtage 
Engine  (Reference  B-12).  The  engine  had  a 040  expansion  ratio  beli  nozzle 
and  delivered  8000  lb^  with  storable  propellants  NjO^/AZ-SO.  Raw  radiance 
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and  transmittance  data  must  be  smoothed  in  order  to  accomplish  the  inversion. 
Deterministic  noise  (instrumentation  fluctuations)  must  be  first  removed  by 
engineering  judgement  and  filtering.  Then,  the  data  are  smoothed  in  a least 
squares  sense.  Figures  B-l  ka  and  B-)  kb  are  raw  and  smoothed  data  for  the 
subject  example.  Figures  B-l  5a  and  B-l  5b  are  inverted  temperature  and  C02 
partial  pressure  radial  profiles. 

Uncertainty 

Random  error  propagation  through  the  inversion  from  the  measured  radiance 
and  transmittance  is  integral  to  one  version  of  the  inversion  program 
(Reference  B- 11).  Typically,  for  a +2  percent,  la  uncertainty  in  radiance 
and  transmittance  in  the  C02  band,  the  rms  uncertainty  in  temperature  is 
1 percent  and  the  uncertainty  in  the  partial  pressure  of  C02  is  6.5  percent. 

Bias  uncertainties  include  calibration  errors  and  unknown  errors  in  band 
model  parameters. 

, 


APPENDIX  B-2.  UV  RESONANCE  ABSORPTION  (REFERENCE  B-12). 


Principle 

If,  in  Figure  B-2  1,  the  transmittance  (t^)  within  a frequency  interval  (Av) 

of  an  electronic  resonance  band  is  known  (measured),  for  each  concentric 

zone  of  constant  temperature  (Tm) , pressure  (Pj  and  species  density  (Nj 

and  T and  P are  known  a priori,  then  the  N are  determinate  from  successive, 
mm  m 

iterate  inversion  of  the  following  set  of  radiative  transfer  equations  using 
the  computational  scheme  of  Figure  B-2  2 (Reference  B-13). 


where 


(tk}Av  " (|r^1  Tm,k  Vk^v 


(B-2  1) 
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(B-2  2) 


j is  incident  spectral  lines; 


i is  absorption  1 ines; 


Av  is  frequency  interval  containing  the  j incident  spectral  lines; 


|°  is  intensity  distribution  of  jth  incident  spectral  line 
Vj 

(assumed  to  be  Doppler  broadened); 


l , is  the  characteristic  length  of  the  m,  k zone; 
m,k 


k is  the  absorption  coefficient  of  the  ith  absorption  line, 

i 

given  by 
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(B-2  3) 


i i a + (u.-y) 


a is  the  ratio  of  coiiisional  and  Doppler  line  half  widths, 
given  by 


a = C P/T  ; C = exp.  constant 
mm 


w.  is  the  Doppler  frequency  function,  given  by 


l\fcn 2 (v.  - v°) 


Wd 


(B-2  1») 


v?  is  the  frequency  at  line  center; 


(^av-)[j  is  the  Doppler  half  width  of  absorber  line  i,  given  by 


<V|>D  ‘ 2 


X is  Boltzmann's  constant; 


c i s velocity  of  1 ight; 


M is  mass  of  absorbing  molecule; 


kyo  is  absorption  coefficient  at  line  center,  given  by 


. 2 Vtt  In  2 e2  , 

v°  T~,  ! Ni"’i 
1 Mc  ^av i D J 


-91- 


/ 


e is  eiectronic  charge; 


Nj,,  is  density  of  absorbing  molecules  in  J"  quantum  state;  for 


equi 1 ibrium, 


Nj„  = No 


9 j. i exp  (-Ejtl/XT) 
fpf) 


N is  total  density  of  absorbing  molecules; 

° I 

; 

g is  statistical  weight  of  J"th  quantum  state; 

J 

Ej,,  is  energy  of  J‘'th  quantum  state; 

(pf)  is  partition  function  for  absorbing  molecules;  and 
f.  is  oscillate  strength  for  the  ith  absorption  line. 

Example 

j 

i 

NO  number  density  profiles  (Figure  B-2  3)  were  obtained  using  absorption  qf 
X-bands  (y=226nm)  (Figure  B-2  4)  through  exhaust  of  turbojet  engine  combustor 
(Reference  B-14).  Temperature  and  pressure  were  measured  separately  by  probes. 
Also  shown  are  results  of  probe  measurements  of  NO  concentration,  indicating 
a well  documented  discrepancy  (References  B-14  and  B- 15). 


Restrictions  and  Weaknesses 


1.  The  technique  is  limited  to  those  molecules  and  atoms  which  have  a 
resonance  transition  in  an  accessible  spectral  range.  The  species 
present  in  combustion  gases  which  have  an  accessible  resonance 
absorption  are: 
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MICROCOPY  RESOLUTION  TEST  CHA 

NATIONAL  BUREAU  OF  ST*E<DAROS-1963- 


TRANSMISSIVITY,  PERCENT  TRANSMISSIVITY,  PERCENT 


Figure  B-2  k.  Measured  Transmittance  at  Peak  of  0,0  -band  (2262 X) 
Through  Exhaust  of  Turbojet  Engine  Combustor  at  Two 
Axial  Locations. 
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766.5,  770  nm 


2.  The  required  inversion  to  obtain  local  values  of  temperature  and/or 
number  density  depends  upon  cylindrical  symmetry  of  the  source. 

3.  The  model  described  above  contains  many  experimental  constants.  The 
technique  is  only  as  good  as  those  constants. 


Implementation 


The  technique  has  been  used  mainly  with  laboratory  equipment,  as  follows 
(Reference  B- 16), 


Source:  specially  constructed  capilliary  discharge  tube 
operating  at  a specified  potential  (3000  vdc) , a 
Specified  pressure  (6  Torr) , and  with  a 12:3:1 


mixture  by  volume  of  A,  N2  and  02»  respectively. 


Receiver:  1/2  M Ebert  spectrometer  operated  at  1.6A  bandpass, 

with  an  RCA  I P28  PM  tube:  data  recording  was  usually 
accomplished  using  strip  chart  type  recorders. 
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A systematic  analysis  of  uncertainty  has  not  been  conducted,  although  an 
approach  similar  to  that  of  Reference  B-ll  Is  quite  adaptable.  A conserva 
tlve  estimate  based  upon  trial  calculation  for  measurement  uncertainty  of 
5 percent  lo  In  the  y-band  of  NO,  the  total  uncertainty  in  number  density 
of  NO  will  be  no  greater  than  + 15  percent  (Reference  B-16). 


1 


APPENDIX  B-3.  RAYLEIGH  SCATTERING  (REFERENCE  B-17). 

Principle 

If,  in  Figure  B-3  I,  the  radiation  spectrum  0V).  scattered  at  an  angle  9 
by  molecules  illuminated  with  a monochromatic  laser  beam  of  intensity  ly 

o 

is  measured,  and  the  relative  species  mole  concentrations  (Cj)  are  known 
a priori,  then  the  gas  density  (N)  and  temperature  (T)  can  be  obtained 
through  the  following  relations  (References  B-17  and  B-18): 


Figure  B-3  1.  Configuration  for  Laser  Rayleigh  Scattering  Measurements 
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where 


K(6)  Is  an  instrumentation  and  orientation  constant; 


o.(v  ) is  the  Rayleigh  scattering  cross  section  for  species  i; 
i o 


ct.  is  the  polar izibil  Ity  of  species  i; 


f . (v)  is  the  velocity  (v)  distribution  function  for  species  i; 


M.  is  the  mass  of  the  Ith  molecular  species; 


X is  Boltzmann's  constant; 


c is  the  velocity  of  light; 


v-v  is  the  Doppler  shift  of  the  scattered  light  due  to  thermal 
o 

motion. 


The  integrated  scattered  intensity  over  the  entire  line.  Equation  (B-3  2), 
is  easily  calibrated  in  terms  of  the  gas  density  if  the  relative  concen- 
trations are  known;  then 


Vv  -CN* 

o 


The  temperature  is  found  by  iterating  on  velocity  until  the  spectral  shape 
of  measured  and  calculated  spectra  (from  Equations  (B-3  l),(B-3  3),  and 
(B-3  *0)  match. 

Restrictions  and  Weaknesses 

I.  Rayleigh  theory  is  valid  only  when  the  scattering  center  is  much,  much 
smaller  than  the  wavelength  of  the  incident  radiation.  Thus,  presence 
of  particulate  matter  destroys  the  technique, and  Hie  scattering 
(Appendix  B- 11)  dominates. 
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2.  The  concentration  (mole  fraction  or  partial  pressure)  of  each  species 
in  the  gas  must  be  known. 

3.  Temperature  measurement  requires  very  high  resolution  (order  of  0.0!  cm  *) 
measurements. 

A.  A density  limit  is  imposed  on  the  temperature  measurement  by  virtue  of 
the  changing  nature  of  the  scattering  from  a Doppler  broadening  to  a 
Brill  ion  broadening,  the  latter  resulting  when  the  molecular  mean  free 
path  becomes  of  the  order  of  the  wavelength  of  the  incident  light 
(Reference  B-  17). 

Implementation 

For  Density  Measurement  (Reference  B-19): 

Laser  source:  I Joule  pulsed  ruby  laser 

Narrow  band  ('v*  IA)  interference  filter  at  laser  wavelength 

Photomultiplier  detector:  standard  electronics 

For  Temperature  Measurement  (Reference  B- 1 7) : 


Laser  source:  argon  ion,  IW  or  more 

-1 

Fabry-Perot  spectrometer:  2 cm  spacing,  0.25  cm  range, 
resolution  0.01*»  cm  ' 


Detector  and  Electronics:  photon  counting  grade  photomultiplier 
tube  with  minicomputer  data  system 


Examples 


1.  Relative  Density  (Reference  B-19).  The  scattering  at  the  incident 
wavelength  from  a pulsed  ruby  laser  beam  passing  through  the  exhaust 
of  500  and  1000  lb^  thrust  model  rocket  engines  burning  storable 
propellants  (N^O^/HHH)  was  measured  at  several  radial  stations.  The 
assumption  was  made  that  the  species  composition  did  not  change  over 
the  radial  extent  of  the  plume;  then  the  relative  gas  density  was 
directly  proportional  to  the  Rayleigh  scattered  intensity.  Knowledge 
of  the  relative  species  densities  was  also  assumed  so  that  the  absolute 
gas  density  profile  was  determined.  Particulate  matter  was  evidently 
very  scarce  in  the  plume  because  the  measured  densities  were  not  very 
different  from  those  calculated  so  that  contributions  from  Mie 
scattering  must  have  been  minimum. 

2.  Temperature  (References  B- 1 7 and  B-20).  The  spectral  line  profile  of 
the  Rayleigh  and  Brill  ion  scattering  from  an  argon  laser  beam  (488  nm) 
passing  through  a hydrogen-air  flame  was  measured  at  several  positions. 
At  the  center  of  the  flame,  temperatures  measured  from  the  line  profile 
matched  well  the  equilibrium  temperature/ and  no  influence  of  Hie  or 
Bril  lion  scattering  was  observed.  At  the  flame  edge  the  Mie  scattering 
from  entrained  particles  from  the  atmosphere  rendered  the  technique 
invalid.  In  the  high  density  precombustion  region  the  Brillion 
scattering  introduced  distortion  of  the  Doppler  broadened  Gaussion 
profile.  The  experiment  provides  opportunity  to  observe  all  the 
phenomena  to  be  encountered  in  the  use  of  Rayleigh  scattering  for 
diagnostics  purposes. 
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APPENDIX  B-4.  RAMAN  SCATTERING  (REFERENCES  B-21  and  B-22). 

Principle 

If,  in  Figure  B-4  I,  the  number  of  photons  within  a wavelength  range  AX 
scattered  per  second  (Q^)  from  a laser  beam  at  wavelength  Xj  traveling 
through  a medium  containing  the  species  8 at  a density  ns  and  a temperature 
T,  due  to  the  Raman  effect  is  measured,  then  ns  (r)  and  T (r)  can  be 
determined  from  the  following  relations: 


RECEIVER 


Figure  B-l*  1.  Illustration  of  Laser  Raman  Scattering  System 


Q*i  * KQ,  ££n  • 

HAX  HX.sj.  sj  sj 


K is  the  instrumentation  constant  which  depends  upon  the  geometry 
collection  efficiency,  and  wavelength  of  the  system; 


Q,  is  the  photon  current  in  the  beam; 


nsj  **  num*)er  of  8 species  molecules  producing  spectral  lines 
j at  wavelengths  within  the  interval  AA  accepted  by  the  receiver; 

°sj  *s  t*,e  Raman  cross  section  for  the  e species  molecules  for 
each  of  the  J spectral  transitions. 


At  thermal  equilibrium 


-e5/xt 

gie  J 


n . - n * i . 
sj  s (pf 


(B-4  2) 


where  g^  is  the  statistical  weight  of  the  initial  energy  level  of  the 
trans i t ion; 


Ej  is  the  energy  of  the  initial  level  of  the  transition; 

X is  Boltzman's  constant; 

j corresponds  to  a given  vibrational  state  v and  rotational  state 
J (for  purely  rotational  transitions  j corresponds  to  a rotational 
state  J) ; 


(pf)s  is  the  partition  function  for  the  s molecule. 

The  Raman  cross  section  for  a specific  vibrational-rotational  transition 
c?n  be  written 


°j  * °v'J'  * Cv'J'  Vv111  Vd1" 


Av'v"J' J" 


where  C is  a constant; 
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gw*  is  the  vibrational  band  strength  factor  relative  to  v'  ««  1, 
v"  * 0 (for  pure  rotational  spectra,  gv'v"  H l); 


Sj , j(l  is  the  relative  rotational  line  strength  factor  (for  diatomic 
molecules,  the  AJ  ■ 0,  ^ 2 Hon  1 -London  factors  are  given  in 
Reference  B-23  for  the  various  angular  momentum  states). 

Ideally,  a Raman  spectrum  from  only  one  molecule  is  radiated  within  AA  so 
that  the  summation  in  s of  Equation  (B-l»  1)  disappears.  Also,  ideally,  only 
one  line  or  band  j of  species  s is  included  in  AA.  Unfortunately,  the  ideal 
is  seldom  met  and  so  the  diagnostic  procedure  contains  a computer  program 
for  determination  of  ng  and  T from  a given  set  of  *AA.  measurement  at 
several  wavelengths  (see  Reference  B-21  for  a detailed  procedure  of  Raman 
Diagnostics) . 

Restrictions  and  Weaknesses 
"3 1 2 

1.  Cross  sections  are  very  small  ('v  10  ^ cm  ) so  that  photon  count  rates 
are  sma 1 1,  and  therefore,  s ignal  to  noise  ratio  (S/N)  is  small.  The 
development  of  the  technique  is  therefore  mainly  an  instrumentation 
problem. 

2.  Because  of  the  inherently  poor  S/N,  spurious  background  radiation 
(typically  from  scattering  off  surfaces,  laser  excited  fluorescence, 
radiation  from  gases  or  particles,  etc.)  limits  applications.  Also, 
the  measurement  time  required  to  achieve  acceptable  S/N  extends  rocket 
firing  times,  increasing  cost  of  testing. 

3.  Overlap  of  lines  from  different  molecules  (s  > I)  within  a practical 
AA  required  to  achieve  a reasonable  S/N  frequently  occurs.  In  rocket 
exhausts  overlap  of  N2  and  CO  rotational  lines  makes  it  necessary  to 
measure  the  combined  density  (Reference  B-23). 
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4.  The  species  that  can  be  measured  are  limited  to  those  gases  having  a 
Raman  transition.  Generally,  this  means  that  atoms  such  as  Na,  K,  0, 
H,  etc.,  are  inaccessible  by  Raman. 

Implementation  (Typical,  Reference  B-23) 

Laser  Source:  Q-switched  pulsed  ruby,  1 Joule  or  more 

Receiver:  Double  tandem  3/4  m spectrometer  with  photon  counting 
quality  photomultiplier  tube  and  minicomputer  data 
system 


Example  (Reference  B-24) 

The  instantaneous  20  nanosec)  temperature  and  species  (N2  + CO  + 02) 
number  densities  on  the  axis  (6  ft  downstream)  of  the  exhaust  plume  from  a 
model  (500  Ib^  thrust),  storable  propellant  rocket  engine  under  simulated 
flight  conditions  was  measured  from  Raman  scattered  radiation  using  a 
pulsed  ruby  laser.  The  purpose  of  the  measurement  was  to  determine  the 
magnitude  of  the  temperature  and  density  fluctuations  in  the  highly  turbu- 
lent mixing  of  plume  and  external  flow.  Figure  B-4  2 is  an  example  of  data 
obtained.  Shown  is  the  rotational  temperature  fluctuations  about  an 
average  value,  the  total  density  fluctuations  of  (N2  + CO  + 02)  about  an 
average  value  as  measured  in  the  rotational  bands  and  the  N2  number  density 
measured  from  the  vibrational  band.  The  AX  for  this  measurement  situation 
was  0.9  nm.  The  firing  rate  of  the  laser  was  12  pulses  per  minute.  The 
uncertainties  quoted  are  based  upon  the  total  number  of  photon  counts 
received. 
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APPENDIX  B-5.  HIGH  RESOLUTION  EMISSION  AND/OR  ABSORPTION 
(REFERENCES  B-25  AND  B-26). 


Principle 

If  individual  spectral  lines  of  a well  defined  vibration-rotation  band  or 
an  electronic  transition  are  resolved,  then  both  temperature  and  number 
density  can  be  determined  from  measured  emission  or  absorption  coefficients. 
The  general  radiative  transfer  equation  through  a path  in  the  gas  stream  is 
written 


•v " exp  ('£\dx)  + exp  (-^Lkvdx)dx 


where  I is  the  monochromatic  radiance  measured; 


1°  is  the  monochromatic  incident  radiance  from  a source; 


k,  is  the  monochromatic  absorption  coefficient  such  that 


/ k d Nf 

o v v 2 m r 
me 


with  e the  electron  charge. 


m the  electron  mass. 


c the  velocity  of  light, 


N the  number  density  of  molecules  in  the  lower 
m 

quantum  state  m. 


f „ the  oscillator  strength  of  the  transition  from 
nm 

the  upper  state  n to  the  lower  state  m; 


(B-5  1) 


(B-5  2) 
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£v  is  the  monochromatic  emission  coefficient  given  by 


dv 


A N hv 
nm  n nm 


(B-5  3) 


with  A the  transition  probability  for  a spontaneous 

nm 

transition  from  state  n to  m, 


N the  number  density  of  molecules  in  the  upper  state  n, 
n 

h the  Planck  constant, 


v the  frequency  of  the  transition  from  upper  state  n 
nm 

to  the  lower  state  m. 


The  oscillator  strength  fnm  is  related  to  the  transition  probability  through 


f 

nm 


it  g N 

me  gn  , . m _n  \ A 

.2  2.2  gm  ' 9n  Nm  nm 


8n  e v_ 


nm 


(B-5  4) 


where  g and  g are  the  statistical  weights  of  state  m and  n respectively, 
m n 

At  equilibrium  temperature  T, 


g -E  /kT 
N =.  m n £ n 
Nn  o „ -Ei/kT 


E g je 


all  i 


g -E  /kT 

M — u ***  ^ 

m o “ .-Ei/kT 


E g.e 


all 


(B-5  5) 


(B-5  6) 
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The  line  shapes  may  be,  generally,  expressed  as  Voight  proviles. 


• y 

V kv  " ev  * kv  7C  ,V..  dy 

nm  nm  a2  + (to-y)2 

nm 


where 


(B-5  7) 


a is  the  ratio  of  the  half  widths  of  lines  due  to  collisional 
nm 

broadening  alone  and  due  to  Doppler  broadening  alone,  and  given  by 


a ■ C £ 
nm  nm  T 


with  C an  experimentally  determined  constant 
nm 


(B-5  8) 


p the  pressure, 


T the  temperature; 


u>  is  the  Doppler  function  given  by 


a = >.■  <n-m-  2y/ln2 


(B-5  9) 


with  ^aV^l/2*  ttie  fre9uer,cy  width  of  the  line  at  half 

intens i ty; 

y is  a dummy  variable  of  integration. 

This  set  of  Equations  (B-5  1)  through  (B-5  9)  is  sufficient  to  describe  any 
spectral  line  transmission  problem. 
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Special  cases: 


Optically  Thin  Case:  kv  “ 0,  I®  * 0 


Equation  (B-5  l)  becomes 


I = jTLe  dx  (B-5  10) 

V O V 

If  the  total  line  intensity  is  measured,  ly  I , S™  dv  = enm, 

in  the  cyl indr ical ly  symmetric  case  e is  obtained  from  the  Abel 

nm 

integral  transform  (Reference  B-27,  for  example) 


e (r)  * - — / 
nm  it  r 


1 -R  dy  *nn/y^ 


(B-5  ID 


/2  2 

/y  -r 


from  e for  two  or  more  lines  one  obtains  local  T and  N . 
nm  o 

Absorption  of  External  Source:  Right  term  in  Equation  (B-5  1) 
eliminated  by  instrumentation  approach.  Then 


'v  " 'C  exp  (-  fo  kvdx) 


a.  1°  * const,  I*;  measure  / -rr  dv  = t , transmittance 
v • nm 

(Reference  B-28) . Then 


ln  T„m  * " (U* 
nm  nm 


(B-5  12) 


(B-5  13) 


For  cylindrical  symetry,  Abel  inversion  again  is  possible,  and 


knm{r)  " • * { 


1 rR  dy  [ln  Tnm(y)] 


(B-5  I1*) 
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Now  local  T and  N„  available  from  k (r). 

o nm 


b.  1°  « 6 (1°  - 1°  ),  v - line  center  frequency  (Reference 

v v v nm 


B-27).  Then 


t(v)  - - /Lk  (v  ) dx 
nm  o nm 


(B-5  15) 


and  again  for  cylindrical  symmetry  local  values  of  k(vnm)  and 

hence  T and  N are  determined  from  (Reference  B-26) 
o 


k(v  ) 
nm 


2* 


\«ta2 


«C 


fnm 

'V'l/2 


(B-5  16) 


This  is  the  "narrow  line  absorption"  case  amenable  to  use  of 
tunable  lasers. 


1°  - 1°  exp 
'*  nm  r 


(B-5  17) 


as  from  a resonance  lamp  where  only  Doppler  broadening  of  the 
source  lines  is  appreciable.  Then 


*Av  " ;Av  {L  ‘v  exP  ( ' l\dx)  dx  dv 


(B-5  18) 


This  is  equivalent  to  the  resonance  line  absorption  technique 
(Reference  B- 1 4)  of  Appendix  B-2  and  the  "onion  peel"  inver- 
sion approach  must  be  taken. 


3.  Line  Reversal,  I®  * Planck  Function, 


Vary  until  ly  ■ (Reference  B-30).  Then 


(Tg)  [1  - exp  (-/kv)dx)  - £Lev(TG)  exP  ("£Lkvdx)  dx  ^B"5  ,9^ 


For  the  constant  temperature,  constant  density  case,  the  integral 
on  the  right  is  the  same  as  the  left  side  and 


Jv<V  ■ cv<TG)L  (8'5  201 


Therefore  TB  ■ T-  and  the  temperature  is  determined.  For  the 

D U 

nonhomogeneous  case  the  equation  must  be  solved  directly  and 
little  advantage  is  served  over  the  general  case,  Equations 
(B-5  1)  through  ( B- 5 9).  Some  sort  of  inversion  must  then  be 
accomplished  in  order  to  obtain  local  properties. 


APPENDIX  B-6.  E-BEAM  FLUORESCENCE  (REFERENCE  B-32). 


Principle 


If  a beam  of  electrons  is  passed  through  a gas,  some  of  the  species  will  be 
col  1 isional ly  excited  to  electronic  states  which  then  radiate  in  some  select 
band  (Figure  B-6  1).  The  radiation  from  the  band  may  then  yield  information 
on  species  density  and  temperature.  At  steady  state,  the  number  of  molecules 
excited  to  a state  n per  unit  volume  per  second,  dQ^/dt,  is  equal  to  the  num- 
ber de-excited  by  radiation  (fluorescence)  to  a state  m plus  those  de-excited 
through  collisions  with  molecular  species, 


Eo  N 
en  o 


I (IA  N - EN  N-0  ,7) 
nm  n . n f nf 
n m f 


(B-6  1) 


where  E is  the  flux  of  electron  flow  through  the  beem; 

o is  the  cross  section  for  excitation  by  electron  collisions; 
en 

Nq  is  the  density  of  the  specie  being  excited; 

A is  the  radiative  transition  probability  from  excited  state 
nm 

n to  a lower  state  m; 

N is  the  steady  state  excited  state  number  density; 
n 

Nf  are  the  densities  of  quenching  molecules; 
onf  are  the  quenching  cross  sections; 
v is  the  mean  velocity  of  the  molecules. 
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ELECTRON 

GUN 


TARGET 


Figure  B-6  1.  Configuration  for  E-Beam  Fluorescence  Measurements 


If  the  quenching  term  can  be  neglected,  Equation  (B-6  I)  simplifies  and  the 
measured  radiation  from  a single  transition  is  just 


-E  /kT 
gen 

I = KA  — c~7ZT 

nm  nm  -t./kT 

?g  i e ' 

i 


from  whence  the  local  temperature  can  be  determined  from  the  distribution. 

From  the  integrated  intensity  of  all  the  lines  in  all  the  bands,  emi tted  the 

local  N is  determined, 
o 

The  above  assumes  that  the  ground  state  rotational  distribution  does 
not  change  during  the  exitation-radiation  process. 

The  technique  can  still  be  used  if  quenching  cannot  be  neglected,  but  the 
Np  v and  must  all  be  known. 


Restrictions  and  Weaknesses 


I.  Quenching  limits  the  technique  to  densities  less  than  about  lO1^  cm 


2.  Only  a few  species  exhibit  e-beam  fluorescence  in  accessible  spectral 
regions.  Applicable  molecules  Include  N2,  H2,  H,  CO,  02>  and  0. 


3.  Interaction  of  low  energy  electrons  and  excited  metastable  species  in 
the  sample  volume  with  excited  states  of  the  species  being  studied 
poses  problems  yet  incompletely  solved  (Reference  B- 33) . 


Equipment 


Source:  15  keV  electron  beam. 


Spectrometer  with  resolution  capable  of  individual  spectral  lines 
*2 


(for  N2,  the  N + 1st  Negative  System  requires  0.2  A*  resolution). 


Appl i cat  ions 


Direct  application  to  the  far  field  of  vacuum  plume  expansions  has  not  been 
made,  although  no  severe  implementation  problems  are  contemplated.  Relevant 
applications  have  been  made  to  simulated  reacting  plumes  (Reference  B-35)  and 
to  mixing  of  flows  in  simulations  of  chemical  laser  mixing  flow  systems 
(Reference  B-34).  The  quenching  and  other  interaction  problems  are  discussed 
in  those  references. 
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APPENDIX  B-7.  RESONANCE  FLUORESCENCE  (REFERENCE  B-36) 


Principle 

If  a beam  of  light  at  energy  hv.  and  intensity  1°  is  incident  on  a gas 

Vfcn 

containing  species  having  a resonance  v^,  (Figure  B-7  0 then  according 
to  the  diagram  of  Figure  B-7  2,  the  rate  of  change  of  the  density  of  the 
excited  state  n is  equal  to  the  number  of  In  photons  absorbed  less  the 
number  lost  by  spontaneous  radiation  to  all  lower  states  m,  plus  those 
lost  by  induced  radiation  to  all  .ower  states  m,  plus  those  lost  during 
molecular  collisions  with  foreign  gas  species,  as  shown  by. 


dN 

n 

dt 


In 


[I  A N + E B N I + E N N,o  ,v] 
nm  n nm  n v . n f nf 

m m nm  f 


(B-7  1) 


At  steady  state  dNn/dt  a 0 and  the  number  density  of  excited  states  is 

V 

* Vffn 

N 


N.  1° 

£ % £n 


Z (A  + B I,  )+ZN,ff,v 
' nm  nm  v , f nf 

m nm  f 


(B-7  2) 


and  the  radiant  emission  coefficient  for  some  n,  k transition  is 


nk 


A . N hv  . 
nk  n nk 


hv  . A . 1°  B0  N. 
nk  nk  v.  £n  £. 

Zn 

Z (A  + B I ) + Q 
nm  nm  v 
m nm 


(B-7  3) 


with  Q substituted  for  the  quenching  term.  Many  special  cases  can  be  considered; 
for  example, 


1.  Q < < (Z  A + B I ),  low  density  case  (Reference  B-37). 
nm  nm  v 
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b i 
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; 


For  any  one  transition  nk 


'nk 


V ‘ok  'v 


fcn 


(B-7  *) 


This  case  is  useful  for  measurement  of  N^  which  can  be  related  to 
temperature  of  the  lower  state  1 if  1 is  one  state  in  a ground 
vibration-rotation  state  distribution,  and  to  the  total  density  of 
the  molecules  for  which  the  resonant  transition  occurs. 


2. 


ZB  1°  > > I A + Q,  saturation  case  (References  B-36  and  B-37). 

m nm  V£m  m nm 

For  a 2-level  system  m=£,  and 


Bn*  ‘In  > > V + Q 


Equation  (B-7  3)  reduces  to 


0 

e .-A  »N  h v „ - A „h  v „ / Jlnx  „„  1 

nl  nil  n n&  nil  nil  (g — j Nil  ^ +-  A 57-s- 

nH  nH_^ 


Bnl'v 


nil 


£n r *„t  hV  ai_  Nt[i  - .=*4^  . , 

9"  Bnt  'v  * J 

nil 


(B-7  5) 


Now  if  the  second  term  in  brackets  is  much  less  than  unity,  the  system  is 

fully  saturated,  and  the  emission  coefficient  is  directly  proportional  to 

the  lower  state  density.  If  the  term  is  not  negligible,  then  a plot  of 

e „ versus  1/1?  will  yield  a straight  line  with  intercept  proportional 
nx  xn 

to  Nj  and  slope  which  yields  the  quenching  rate,  Q. 


The  measured  intensity  of  the  radiation  scattered  into  some  detector  is 
proportional  to  the  emission  coefficient, 
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Restrictions  and  Weaknesses 


1.  The  technique  can  be  used  in  a lower  and  upper  density  range,  but  not 

in  between  because  of  the  unknown  quenching  rate.  The  estimated 
1 6 —3  19  “3 

limits  are  N < 10  cm  J or  N > 10  7cm  3 corresponding  to  negligible 
and  saturation  quenching  rates,  respectively. 

2.  Limited  to  molecules  or  atoms  having  an  electronic  resonance  transi- 
tion: OH,  CN,  CH,  C^,  Na,  K,  etc. 


Equipment 

Source:  tunable  dye  laser  in  the  resonance  wavelength  range 
of  the  species  of  interest. 

Receiver:  spectrometer  or  interference  filter  with  photomultiplier 
detection  system. 

Examp  1 e 

The  only  reported  example  akin  to  rocket  flows  is  the  concentration 
measurement  by  saturated  fluorescence  in  an  atomospher ic  pressure  oxygen/ 
acetylene  flame  reported  in  Reference  B-37.  The  value  found  was  4.5X10^ 
cm  ^ with  an  estimated  uncertainty  of  300  percent.  The  large  uncertainty 
resulted  mostly  from  lack  of  accurate  calibration  in  the  proof  of  principle 
experiment  and  is  not  considered  representative  of  the  technique. 
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APPENDIX  B-8.  HOLOGRAPHY  (REFERENCE  B-38) 


Principle 


If  the  undisturbed  amplitude  of  a monochromatic,  coherent  electromagnetic 
field  is  mixed  with  the  amplitude  of  a wave  diffracted  by  an  object  in  the 
same  coherent  field,  the  result  will  be  a diffraction  pattern,  or  inter- 
ferogram,  containing  the  Fourier  transform  of  the  disturbed  wave.  The 
Fourier  transform,  which  contains  three-dimensional  information  about  the 
object,  can  be  recorded  on  a photographic  emulsion.  If  the  emulsion  is 
exposed  in  a linear  portion  of  the  film  density  versus  intensity  relation, 
then  the  film  density  represents  an  accurate  recording  of  the  Fourier 
transform.  The  inverse  of  the  Fourier  transform  again  represents  the 
object  field  and  is  accomplished  by  exposing  the  interferogram  to  mono- 
chrometric,  coherent  light  and  reforming  the  image  of  the  object.  This 
process  is  referred  to  as  holography.  The  literature  describing  the  many 
aspects  of  holography  is  vast.  Reference  B- 38  describes  the  use  by  holo- 
graphy in  flow  visualization  and  particle  measurement  application. 


If  the  objective  field  is  a particle  cloud  then  an  image  of  the  reconstructed 
hologram  can  provide  a magnified  cross  section  of  the  particle  field, 
yielding  size  and  shape  of  each  particle  that  is  large  enough  to  be  imaged 
by  the  optical  system.  If  a double  pulsed  laser  is  used,  the  movement 
of  a particle  between  pulses  can  be  measured  by  a double  exposure  hologram 
and  the  velocity  obtained.  The  change  in  size  and  shape  between  pulses 
yields  information  on  such  things  as  growth  rate  in  a condensing  flow,  or 
surface  reaction  rate  during  combustion.  If  the  particle  is  moving  very 
fast,  then  the  film  streak  formed  during  a pulse  interval  (or  a pulse  chain) 
of  the  laser  provides  a measure  of  both  velocity  and  size,  although  some 
degradation  in  resolution  is  suffered. 


If  the  objective  field  is  a variation  in  index  of  refraction,  or  density, 
then  the  illuminating  beam  may  be  mixed  with  a reference  beam  in  an  inter- 
ferometer arrangement.  The  resulting  interferogram  may  then  be  recorded 
as  a hologram  and  analyzed  in  terms  of  flow  field  density  variation. 


1 

i.  I 

1 1 
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Restrictions  and  Weaknesses 


1.  The  particle  size  is  limited  by  the  ordinary  diffraction  limit  of  all 
optical  systems  to  a smallest  dimension  R 

R - 1.3  A f 

where  X is  the  wavelength  of  the  incident  light,  and  f is  the  ratio  of 
the  usable  cone  of  light  at  the  optical  element  (lens)  to  the  distance 
from  the  object  to  the  plane  of  the  optical  element.  For  a ruby  laser 
and  an  f2  system,  the  best  resolution  attainable  would  be  about  2 ym. 
The  depth  of  field  Z,  over  which  the  hologram  is  approximately  in 
focus  given  by  (Reference  6*38) 

d2 

Z - .00  f 

represents  another  limit  such  that  the  field  cannot  be  defined  well 
outside  that  interval.  For  a ruby  laser  and  a d * 2 ym  particle, 
the  depth  of  usable  field  would  be  about  600  ym.  Imperfections  in 
optics,  noise  from  diffraction  due  to  gas  density  gradients,  etc., 
usually  prevent  these  theoretical  limits  from  being  obtained. 

2.  Data  reduction  is  a major  problem  because  of  the  difficulty  of  placing 
quantitative  criteria  on  focus  quality  and  noise  suitable  for  mechaniza 
tion.  Consequently,  human  judgement  is  generally  required  so  that  data 
reduction  of  a flow  field  with  heavy  particle  loading  is  time  consuming 

3.  "Freezing"  of  particles  require  very  short  pulse  lasers.  In  order  to 
prevent  blurring  of  particles  the  pulse  width  should  be  less  than 
about  0.1  d/v  where  d is  the  particle  diameter  and  v is  the  particle 
velocity.  Thus  for  a 2y  particle  traveling  at  only  200  m/sec,  the 
pulse  width  should  be  less  than  I nanosec.  Commercial  lasers  are 
currently  limited  to  about  20  nanosec,  thus  limiting  velocities  at 
which  2 ym  resolution  holography  can  be  applied  to  about  20  m/sec. 


For  this  reason,  holography  appears  to  be  limited  to  combustion  chamber 
phenomena  or  the  far  field  condensation  problem. 

A.  Velocity  measurement  by  double  pulse  holography  is  limited  by  the  time 
between  laser  pulses  relative  to  the  particle  size  and  the  fact  that 
the  particle  must  be  distinguishable  at  the  two  locations.  For  a 10  pm 
particle,  a rule  of  thumb  that  the  particle  should  not  move  more  than 
about  10  diameters  and  a 0.5  psec  pulse  interval  (state  of  the  art  of 
commercial  lasers).  The  maximum  velocity  that  can  be  measured  is  about 
200  m/sec  (Reference  B-38).  This  criteria  is  quite  subjective  but 
suggests  that  velocities  no  larger  than  about  200  m/sec  are  practical 
for  particles  less  than  10  pm  dia.  Velocity  measurement  from  photo- 
graphic streaks  or  multiple  exposures  are  limited  basically  by  the 
depth  of  field  of  the  optical  system.  For  the  600  pm  field  and  a 2 
nanosec  pulse  the  velocity  would  be  30  km/sec,  but  film  exposure  limits 
prevent  this  limit  from  being  obtained. 

Implementation 

A great  variety  of  optical  configurations  have  been  used  for  obtaining  and 
reconstructing  holograms.  For  applications  in  combustion  chambers  the 
system  described  in  Reference  B-39  and  shown  schematically  in  Figure  B-8  1 
seems  near  optimum. 


Example  (Reference  B-39) 

The  system  of  Figure  B-8  1 was  used  to  obtain  holograms  from  a windowed 
combustion  bomb  containing  different  types  of  solid  propellants.  The  2- 
inch  interval  diameter  bomb  was  equipped  with  opposing  3/4-inch  thick 
windows  and  was  operated  at  either  500  or  1000  psi  chamber  pressure  with 
a firing  duration  of  about  300  millisec.  The  laser  had  a pulse  width  of 
about  50  to  80  nanosec  with  a separation  between  pulses  of  about  1.3  psec. 
With  the  diffuser  plate  in  Figure  B-8  1 inserted,  good  quality  holograms 
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LENS-PLATE  STRUCTURE  USED 


Arrangement  for  Application  of  Holography  to 
issure  Solid  Propellant  Combustion  Studies 


1 


were  recorded  with  resolution  down  to  about  2.5  Mm.  With  the  diffuser 
plate  removed,  the  photographs  showed  only  the  refraction  patterns  due  to 
changing  indices  of  refraction  in  the  gaseous  combustion  zones.  A 
particle  size  distribution  obtained  with  this  system  is  shown  in  Figure 
B-8  2.  The  data  reduction  was  accomplished  manually.  Velocities  on 
the  order  of  10  m/sec  were  measured. 


Figure  B-8  2.  Example  Particle  Size  Distribution  Data  from  Combustion  Bomb 

Using  Solid  Propellant  Sample  at  500  psia.  Obtained  with  Holo- 
system  of  Figure  B-8  1 (Reference  B'39). 


APPENDIX  B-9.  LASER  DOPPLER  VELOCIMETRY  (LDV) (REFERENCE  B-38) 

Principle 

The  Mie  scattered  light  from  a coherent,  monochromatic  beam  (laser)  incident 
on  a small  particle  moving  at  a velocity  v will  be  shifted  by  frequency  by 

AVD  " VD  ' Vo  " Hf“  cos6  (B"9  0 

where  vQ  is  the  Doppler  frequency 

v Is  the  particle  velocity 

A is  the  laser  wavelength 

9 is  the  angle  between  the  observer  and  the  particle  velocity  vector. 

Thus,  measurement  of  A provides  a measure  of  the  velocity.  Modern  Doppler 
velocimetry  (or  anemometry)  revolves  around  methods  to  measure  A Vp. 

1.  Direct  Measurement:  The  shift  in  the  frequency  of  the  light  beam  is 
measured  directly, 

v.  * v + — r—  cos0  (B-9  2) 

D O A 

using  a Fabry-Perot  interferometer  (Reference  B-Ao).  (See  Figure  B-9  1.) 

2.  Reference  Beam:  Some  part  of  the  incident  beam  is  split  off  and  mixed 
with  the  scattered  beam  and  added  on  a detector.  The  result  is  an 
optical  heterodyning  signal,  I,  such  that, 

I * Aj  + A^  + 2Aj  A2  cos  2t r (vD  - v^Jt  (B-9  3) 

where  Aj  and  A2  are  the  amplitudes  of  the  two  scattered  signals  and 
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is  the  frequency  difference,  A v^,  the  Doppler  frequency.  The  frequency 
is  easily  determined  using  conventional  electronic  instrumentation  (Ref- 
erence B-41)  (see  Figure  B-9  2). 

3.  Dual  Scatter:  Two  beams  are  incident  on  the  scatterer  with  an  angle  <p 
between  the  beams,  and  the  scattering  from  both  beams  is  summed  on  a 
detector  mounted  at  any  orientation  with  the  scattering  volume.  See  Figure 
B-9  3.  The  difference  heterodyning  frequency  is  again  Avp  , now  related 
to  the  velocity  through 

A VD  " 2 tT  ^1  ' V"*  ^ 

which  can  be  written, 

A vD  - ~ sin  (B-9  5) 

as  shown  in  Figure  B-9  k (Reference  B-42). 


I*.  Fringe  Model  (Velocity):  In  the  dual  scatter  configuration  (Figure 

B-9  5),  the  two  interesting  beams  actually  form  an  interference  fringe 
pattern  in  the  focal  volume  with  spacing 


. A 

2s in  <t> 
2 


(B-9  6) 


Hence,  a particle  with  diameter  d less  than  6,  in  crossing  the  fringe 
pattern  at  velocity  v,  produces  an  oscillating  scattering  signal  with 
frequency  f equal  to, 


(8-9  7) 


Now  f is  Just  A vD  of  the  dual  scatter  method,  and  thus  the  two  inter- 
pretations of  the  configuration  are  equivalent  (Reference  B-l»3). 


5.  Fringe  Model  (Size):  The  shape  of  a scattering  event  as  the  particle  of 
diameter  d traverses  the  fringe  pattern  of  spacing  6 is  related  to  the 
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particle  size  through  the  "Visibility",  defined  by  (Reference  B-1*1*) 


v 


max  m i n 

I + I . 
max  m i n 


2 J,(x) 
x 


where  X is  d/6  and  Jj  is  a Bessel's  function  of  the  first  kind.  In 
practice,  a calibration  is  required  since  the  proportionality  constant 
is  not  well  known  and  depends  upon  the  optical  configuration  (Ref- 
erence B-l»5.) 


Restrictions  and  Weaknesses 

1.  An  upper  limit  on  velocity  using  electronic  processing  is  set  by  a 

g 

practical  upper  frequency  limit  on  electronic  gear  of  about  10  Hz. 

This  translates  into  upper  velocity  limits  of  a few  hundred  meters 
per  second  for  the  direct  heterodying  method,  depending  upon  the 
angle  of  observation,  and  about  1.5  km/sec  for  the  dual  scatter 
(fringe)  technique,  depending  upon  the  lower  practical  value  of  the 
angle  between  the  intersecting  beams  (assumed  2°  here) (Reference  B-  **6) . 

2.  Fabry-Perot  interferometers  permit  measurements  of  frequency  shifts 
corresponding  to  very  high  velocities,  but  a 10  MHz  resolution  limit 
translates  into  a resolution  uncertainty  in  velocity  of  about  10  m/sec. 
The  greatest  weakness  of  this  method  of  measuring  the  Doppler  shift, 
however,  is  in  the  difficulty  of  maintaining  optical  stability  under 
field  conditions  (Reference  B- 39) . 


3.  LDV  measures  particle  velocity;  the  desired  measurement  is  usually 
gas  velocity.  Viscous  drag  may  affect  particle  velocities,  depending 
upon  gas  density  and  particle  size,  to  an  extent  where  particle 
velocity  lags  gas  velocity.  Thus,  two  phase  flow  theory  must  be 
applied  to  determine  if  this  problem  exists  in  any  specific  applica- 
tion (Reference  B-  ^7) . 
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4.  Measurements  in  turbulent  flow  take  on  their  own  characters  and  must 
be  analyzed  in  terms  of  the  statistical  behavior  of  both  the  particle 
sampling  and  the  turbulent  structure  (Reference  B-47). 

5.  The  range  of  part.cle  sizes  measurable  by  the  shape  of  the  fringe 
pattern  is  limited  on  the  lower  end  by  the  asymptotic  behavior  of  the 
visibility  at  about  d/6  ■ 0.2  (see  Figure  B-9  6);  that  is,  d = 0. 1 A/ 
sin<|>/2,  which  for  <f>  * about  12°  as  a maximum  gives  a minimum  value  of 
about  d = A (Reference  B-48) . The  maximum  value  is  set  at  about  the 
value  of  d/6  at  which  the  visibility  becomes  double  valued  (see  Figure 
B-9  6);  for  spherical  particles  this  limit  is  about  d/6  = 0.8,  which 
for  the  4>  = 12°  system  with  an  argon  laser  places  an  upper  limit  of 
about  2.5  pm  on  the  size.  Thus,  for  a fixed  system  the  range  of  sizes 
measurable  is  not  large  (i.e.,  0.5  < d < 2.5  pm  for  a 12°  argon  laser 
system;  3 < d < 12  pm  for  a 1°  system). 

6.  Interpretation  of  the  scattering  burst  pattern  is  required  to  determine 
whether  the  particle  is  traversing  the  fringe  pattern  perpendicular  to 
the  fringe  pattern  or  at  an  angle,  to  determine  if  more  than  one  particle 
is  in  the  scattering  volume  at  a time,  and  to  determine  any  number  of 
other  eventualities.  This  places  very  stringent  requirements  on  the 
signal  processor,  and  the  data  reduction  and  weakens  the  technique 
(Reference  B-48) . 


Implementation 

Equipment  used  in  LDV  systems  takes  on  many  forms.  A particularly  successful 
system  is  described  in  Reference  B-48 , whi-h  was  used  for  both  velocity  and 
particle  size  measurements  inside  a diesel  engine  combustion  chamber.  A 
system  designed  for  high  velocity  measurements  in  turbojet  engine  exhausts 
is  described  in  Reference  B-49. 


3 

■i 
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Measurement  of  the  axial  velocity  distribution,  the  rms  turbulent  fluctuations 
and  the  turbulent  intensity  made  in  the  exhaust  of  an  F 1 01 -GE- 1 00  turbofan 
engine  are  described  in  Reference  B-^9.  The  dual  scatter  system  utilized 
the  fringe  crossing  principal  and  analyzed  individual  scattering  events  to 
obtain  the  measurements.  Example  data  are  shown  in  Figure  B-9  7.  The  sys- 
tem had  a maximum  velocity  limit  of  about  *»500  ft/sec  (1370  m/sec). 
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ITY  DISTRIBUTION  PROBABILITY  DISTRIBUTION  PROBABILITY  DISTRIBUTION 

ON,  P(V)  FUNCTION,  P(V)  FUNCTION,  P(V) 


MEAN  VELOCITY,  V - 910.9  fps 
rms  VELOCITY,  V'  - 37.5  fps 
TURBULENCE  INTENSITY,  a - 0.041 

o CENTERBAND  P(V) , AV  - 28.9  fps 


) 1,000 

AXIAL  VELOCITY,  V,  fps 


1,500 


MEAN  VELOCITY,  V * 2360.7  fps 
rms  VELOCITY,  V'  104.0  fps 
TURBULENCE  INTENSITY,  c - 0.044 

o CENTERBAND  P(V) , AV  * 45.9 


?,500 


1 \ 


2,000 


2,500 


INTERMEDIATE  POWER 


3,000 


AXIAL  VELOCITY,  V,  fps 


MEAN  VELOCITY,  V * 2973.9  fps 
rms  VELOCITY,  V'  - 201.8  fps 
TURBULENCE  INTENSITY,  o * 0.068 

o CENTERBAND  p(v) , AV  ■ 45.9  fps 


MAX  THRUST  AUGMENTATION 


AXIAL  VELOCITY,  V,  fps 

Figure  B-9  7.  Velocity  Histograms  of  the  Axial  Velocity  at 
Nozzle  Exit  of  the  Exhaust  of  an  FI01-GE-I00 
Turbofan  Engine  Obtained  by  Laser  Doppler  Velocimetry 


APPENDIX  B-10.  DIFFUSION  CORRELATED  SPECTROSCOPY  (DCS) 
(REFERENCES  B-50,  B-51,  B-52) 


r 


Principle 


If  the  power  spectrum  at  low  frequency  (id  ■ 0 to  ~20  kHz)  of  the  light 
scattered  by  monodisperse  particles  In  Brownian  motion  at  an  angle  4> 
from  a uniform  illumination  source  (laser)  is  measured  (Figure  B-10  1), 
spectrum  of  the  Lorentzian  form 


Ps(u>) 


K2D/tt 

2 2 2 

,KD)V 


(B-10  1) 


will  be  observed  (Figure  B-10  2), 

where  K is  the  absolute  value  of  the  scattering  vector 


K - -S-  sin  on-,  I®'10  2> 

o 

\ is  the  wavelength  of  the  incident  laser  beam; 
o 

D is  the  diffusion  coefficient  of  the  particles,  which 
for  Stokes  flow  (molecular  mean  free  path  small  with 
respect  to  particle  size)  is  given  by 


n - k T c (B-10  3) 

° 3irnd 

k is  Boltzmann's  constant; 

T is  the  temperature; 
n is  the  coefficient  of  viscosity; 
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C is  a slip  flow  correction  factor  (for  a molecular  mean  free 
path  small  compared  to  the  particle  size,  C -*•  1); 

d is  the  particle  diameter. 


The  frequency  in  Hz  at  the  half  power  point  is 


W 


1/2 


2 JL. 

2 


sin 


Jl__L 

2 d 


(B-10  A) 


so  that  either  T or  d may  be  determined  from  the  measured  half  width. 


This  result  makes  a powerful  tool  of  a rather  simple  measurement  scenario. 
The  result  is  based  upon  an  involved  theoretical  analysis  of  the  scattering 
from  randomly  spaced  particles  in  Brownian  motion  (given  completely  in  Ref- 
erence B-53  and  in  considerable  detail  in  Reference  B-53  and  8-5*0-  If  the 
illumination  field  is  not  uniform,  then  the  intensity  distribution  along 
the  path  of  observation  must  be  convoluted  with  the  phase  function  produced 
by  the  Brownian  motion,  and  a power  spectrum  of  Voight  form  results  when 
the  distribution  is  Gaussian.  The  latter  should  be  avoided  by  proper 
experimental  arrangement,  if  possible,  and  further  consideration  of  that 
case  (treated  in  References  B-51  and  B-5*»)  will  not  be  given  here. 


Restrictions  and  Weaknesses 

1.  The  technique  is  considered  to  be  still  in  the  developmental  stage. 
Two  important  extensions  that  must  be  made  are  (1)  development  of  the 
mathematical  model  to  invert  a set  of  measurements  made  at  different 
angular  orientation  to  the  parameters  of  a polydisperse  distribution 
of  sizes,  and  (2)  introduction  of  a more  adequate  particle  drag  law 
to  handle  lower  density  flows  where  the  Stokes  law  is  not  applicable. 


2.  A fundamental  limitation  on  the  particle  size  range  measurable  at  any 
given  velocity  exists  because  of  a limitation  on  observation  time. 

One  way  of  looking  at  this  limitation  is  to  say  that  sufficient 
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observation  time  must  exist  for  the  random  Brownian  motion  to  assert 
itself;  in  a similar  manner,  a time  interval  (or  a sufficient  number 
of  observations)  is  required  in  the  classical  oil  drop  experiment  to 
establish  the  Gaussian  distribution  for  Brownian  motion.  Penner, 
Bernard,  and  Chany  (Reference  B-55)  have  considered  this  problem  and 
have  determined  the  effect  of  including  a finite  observation  time  T. 
Their  result  can  be  stated:  the  Lorentzian  spectral  function  will 
not  be  perturbed  appreciably  if. 


>0.2  (-^-) 
T 


or,  where  T * Az/v,  where  Az  is  the  diameter  of  the  laser  beam,  this  criteria 
can  be  written 


vd  < 0.9  Az 

Ap 

with  u the  viscosity,  k Boltzmann's  constant,  and  T the  temperature.  For  a 
flame  at  2000°  K,  a He-Ne  laser  (A  ■ 0.6328  ym)  and  y ■ 8 X 10  ^ poise, 

vd  < 0.08Az 

A few  representative  values  follow  for  illustration. 


d(ym) 

v (m/sec) 

Az(cm) 

0.1 

1 

0.0125 

0.1 

10 

.125 

0.1 

100 

1.25 

0.1 

1000 

12.5 

1 

10 

1.25 

.01 

100 

0.125 

Thus,  as  velocity  increases  the  illumination  volume  must  be  increased,  which 
introduces  spatial  resolution  limitations  and  requires  larger  power  lasers 
to  retain  signal  to  noise  ratio. 
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Implementation 


Figure  B-10  3 shows  the  essential  equipment  required  for  the  diffusion 
correlated  spectroscopy  technique.  A typical  experiment  conducted  at  the 
AFRPL  (Reference  B-56)  utilized  the  following  equipment  : 75  mw  He-Ne 

O © 

laser,  1.6A  bandwidth  centered  at  6328A,  RCA  1*832  PH  tube,  tape  recorder 
(bandwidth  -25  KHz  fm  systems).  Spectral  Dynamics  spectrum  analyzer 
Model  360-65. 


Example  (Reference  B-56) 

The  system  shown  in  Figure  B-10  1*  was  used  to  obtain  data  in  an  acetylene/ 
air  flat  flame  burner.  A spectrum  is  shown  in  Figure  B-10  k with  the  data 
reduction  indicated,  giving  a particle  size  of  about  312  nm. 


Apparatus  for  Laser  Doppler  Line  Shape  Measurement. 


(Scattering  Angle  - 60°,  Height  Above  Burner  Face  - 3 cm,  Equivalence  Ratio  - 2.5) 


^ 


fi 

I 
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APPENDIX  B-H.  MIE  SCATTERING  INVERSION  (REFERENCE  B-57) 

Principle 

The  steady  state  angular  distribution  of  scattered  power  from  a laser  beam 
passing  through  a medium  containing  particles  (for  single  scattering)  (Figure 
8-11  1)  is  given  by  (Reference  B-57) 

P(6)  * C Ja2  [i y (e,m,a)  + i2  (0  ,m,«)J  f (a)  da  (B-11  1) 

ai 

where  P(0)  is  the  scattered  power; 

6 is  the  angle  of  observation  with  respect  to  the  laser  axis; 

C is  a constant  dependent  only  upon  the  experimental  configuration; 

ij  is  the  Mie  scattering  function  for  parallel  polarization  (defined 
in  Reference  B-58,  page  47); 

i2  is  the  Hie  scattering  function  for  perpendicular  polarization, 
(defined  in  Reference  B-58,  page  47); 

m is  the  index  of  refraction  of  the  particles; 

a(*nd/X)  is  the  size  parameter;  d is  the  diameter  of  the  particles 
and  X the  laser  wavelength; 

f(a)  is  the  particle  size  distribution  function; 
oij  ,a2  define  the  range  of  a. 
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P^  (cos  0)  are  associated  Legendre  Polynomials  of  the  first  kind; 
An(m,o)  and  bn(m,a)  are  linear  combinations  of  Ricatt i-Bessei 
functions. 
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For  a given  set  of  6,  m and  a the  scattering  functions  can  be  calculated 
in  a straightforward  manner  by  a high  speed  computer.  The  diagnostics 
problem  then  becomes  a matter  of  inverting  the  integral  Equation  (B-l I I) 
to  determine  a distribution  of  a *rom  the  measured  P(6). 

The  problem  can  be  cast  in  a more  manageable  form  by  choosing  a finite 
number  of  measurement  angles  0j,  and  a finite  number  of  particle  size 
intervals  Aa^,  then  in  matrix  form  [with  P(0.)/C  = G (6 . ) ] 

G - SF  (B-ll  6) 

where  G is  a column  vector  of  the  measurements  at  the  6.,  F is  a column 
vector  of  the  number  of  particles  within  each  size  element  Aa^,  and  S is 
the  square  transformation  matrix  with  elements  e . ^(0 . ,m,ak)  calculable 
from  the  scattering  functions  if  the  index  of  refraction  is  assumed.  The 
inverse  of  the  matrix  Equation  (B-11  6)  constitutes  the  determination  of 
the  distribution  function  f(otk).  This  inversion  is  far  from  a trivial 
numerical  problem,  as  discussed  in  References  B-57»  B-59,  and  B-60. 

The  amount  of  information  from  a set  of  measurements  may  be  doubled  by 
selectively  measuring  the  scattered  parallel  and  perpendicularly  polarized 
component  of  the  intensity.  Also,  an  additional  piece  of  information  is 
given  by  the  extinction  measurement  at  0*0.  Various  investigators  (References 
d-  Go,  B-61,  and  B-62,  for  example)  have  used  combinations  of  angular, 
polarization,  and  extinction  measurements  to  infer  approximate  particle 
size  distributions. 

I 

Status 

This  technique  is  in  the  research  stage  with  no  successful  general  application 
reported  to  date.  Lewis,  et  al,  at  AEDC  (Reference  B-60)  have  an  operating 
computer  program  which  they  have  used  to  generate  particle  size  histograms 
when  only  two  angular  measurements  could  be  made. 
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Restrictions  and  Weaknesses 


1.  The  complex  nature  of  the  deconvolution  (or  inversion)  process  makes 
this  technique  difficult  to  access  until  some  definitive  application 
is  made.  Even  more  difficult  may  be  the  assignment  of  uncertainty 
limits  to  the  data.  The  uncertainty  problem  is  discussed  in  Reference 
B-57  in  the  context  of  its  contributions  to  instability  of  the  inversion 
calculation.  Problems  of  uniqueness  of  the  inversion  have  also  been 
posed  but  are  as  yet  unresolved. 

2.  The  complex  index  of  refraction  must  be  known  or  must  be  determined 
as  a third  order  iteration  on  the  inversion.  Currently,  very  little 
data  on  the  index  of  refraction  of  materials  of  interest  in  rocket 
exhausts  exist.  This  lack  of  data  on  the  index  of  refraction  may  be 
the  weakest  link  in  the  use  of  the  technique. 

3.  Limits  on  the  size  range  of  applicability  of  the  technique  occur  when 
the  change  in  scattering  function  with  size  becomes  asymptotically 
small.  On  the  lower  end  the  scattering  function  becomes  relatively 
insensitive  to  size  changes  at  a diameter  of  about  0.25A  (Rayleigh 
limit).  On  the  upper  end  the  limit  on  measurable  particle  size  is 
probably  about  10A,  although  this  limit  is  really  a function  of  how 
well  the  back-scattered  radiation  can  be  measured. 

Implementation 

The  apparatus  necessary  for  this  measurement  consists  of  a laser  having 
sufficient  power  to  penetrate  the  particle  laden  flow  and  the  appropriate 
detector  systems.  No  suitable  apparatus  for  rocket  exhausts  has  been 
tested  to  date,  although  a system  is  under  development  at  the  AFRPL. 
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